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Abstract 

This  research  presents  the  analysis  and  verification  of  a  stripline  designed  by 
Air  Force  Research  Laboratory  for  use  in  measuring  the  electrical  properties  of  mate¬ 
rials  at  low  frequencies  and  high  temperature,  ft  is  designed  to  operate  in  the  TEM 
mode  up  to  4  Ghz  and  have  a  characteristic  impedance  of  50  ohms.  A  full  wave 
base  method  is  used  to  analyze  the  structure.  The  parallel  plate  waveguide  dyadic 
Green’s  function  is  developed  for  a  current  immersed  in  a  PEC  parallel  plate  environ¬ 
ment.  ft  is  used  to  formulate  a  pair  of  coupled  electric  field  integral  equations  (C1E). 
These  CIEs  are  solved  through  a  computationally  efficient  entire-domain  method  of 
moments  (MoM)  technique.  Numerical  efficiency  is  gained  through  employing  Cheb- 
syshev  polynomials  of  the  first  and  second  kind  as  testing  and  expansion  functions. 
These  efficient  expansion  and  testing  function  sets  require  as  few  as  three  expansions 
for  accurate  results.  Further  numerical  efficiencies  are  gained  by  taking  advantage 
of  transverse  electromagnetic  propagation  properties  to  develop  a  specialized  TEM 
integral  equation  reducing  the  number  of  integrations  performed  from  2N2  +  N  to  N 2 . 
An  expression  for  the  characteristic  impedance  is  developed  using  the  MoM  results. 
The  characteristic  impedance  is  calculated  for  various  degrees  of  center  conductor 
miss  alignment. 
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Full- Wave  Based  Validation 
of  Stripline  Field  Applicator 
for  Low  Frequency 
Material  Measurements 

I.  Introduction 

1.1  Problem  Statement 

The  objectives  of  this  research  are  to  validate  the  design  parameters  of  a  high 
temperature/low  frequency  stripline  apparatus  developed  for  the  Air  Force  Research 
Laboratory  using  a  full-wave  analysis  and  to  predict  the  effects  of  center  conduc¬ 
tor  misalignment  on  the  characteristic  impedance  of  the  apparatus  used  in  material 
characterization  measurements. 

1.2  Motivation  for  Stripline  Use 

Current  material  testing  at  low  frequencies  requires  waveguides  having  large 
cross-sectional  dimensions  (e.g.  four  feet  wide  by  two  feet  high).  Fabricating  sam¬ 
ples  this  large  is  costly  and  exceedingly  difficult.  In  addition,  obtaining  an  oven  to 
heat  this  waveguide/samplc  fixture  would  require  enormous  space  and  power,  mak¬ 
ing  it  virtually  impossible  to  implement.  The  proposed  stripline  design  would  have 
cross-sectional  dimensions  of  approximately  five  inches  by  two  inches,  thus  samples 
are  easily  fabricated  and  pre-existing  furnaces  can  easily  be  used.  Thus,  the  stripline 
fixture  will,  for  the  first  time,  give  Air  Force  Research  Laboratory  the  capability  of 
measuring  shielding  material  at  low  frequencies  and  high  temperatures.  AFRL’s  ap¬ 
proximate  design  specifications  stipulate  that  a  50  ohm  stripline  can  be  realized  if 
ground  plane  spacing  is  approximately  3.5  cm  and  center  conductor  width  is  approxi¬ 
mately  5  cm.  A  full- wave  analysis  will  be  developed  in  this  research  to  more  precisely 
validate  the  design. 
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1.3  Motivation  for  Investigating  Center  Conductor  Misalignment 

Striplines  used  to  characterize  the  properties  of  materials  are  air-filled  to  al¬ 
low  sample  insertion.  The  unsupported  center-conductor  often  sags  (see  figure  1.1) 
leading  to  characteristic  impedance,  Zq,  uncertainty.  It  has  been  found  that  permit¬ 
tivity  and  permeability  measurements  are  highly  susceptible  to  S-parameter  and  Z0 
uncertainty  [2,18].  Better  characterization  can  be  achieved  with  an  expression  for  Zq 
that  can  accommodate  center-conductor  misalignment.  Previous  techniques  include 
conformal  mapping,  finite  difference  equations  and  various  Green’s  function/integral 
equation  methods.  While  conformal  mapping  yields  exact  results,  it  involves  compli¬ 
cated  functions  and  is  cumbersome  to  implement  especially  when  the  system  is  not 
symmetric.  The  finite  difference  technique  requires  evaluation  points  on  the  order  of 
5,000-10,000,  making  it  computationally  intensive  [25].  Integral  equation  approaches, 
like  Kammler’s,  lower  computational  size  by  using  a  sub-domain  method  of  moments 
technique.  Approximate  methods  are  also  presented  in  the  literature  [20,31].  The 
goals  of  this  research  is  to  develop  a  full-wave,  computationally  efficient  expression 
for  computing  Z0  of  a  stripline  having  a  misaligned  center  conductor. 

1.4  Limitations 

In  order  to  make  the  problem  tractable  a  simplified  model  of  the  stripline  has 
been  assumed.  The  geometry  is  considered  to  be  invariant  along  the  guiding  axis 
and  the  center  conductor  is  modeled  as  an  infinitesimally  thin  strip.  In  addition, 
only  pre-existing  numerical  algorithms  in  Matlab  are  used  for  numerical  integration. 
The  numerical  integration  techniques  implemented  in  this  research  are  an  adaptive 
Simpson  quadrature,  the  trapezoidal  method,  an  adaptive  Lobatto  quadrature,  and 
adaptive  Gauss-Kronrod  quadrature.  These  various  methods  are  used  to  accommo¬ 
date  the  various  forms  of  integrals  encountered  in  this  research. 
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Figure  1.1:  Example  of  center  conductor  sag. 
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1 . 5  Scope 

There  are  certainly  other  methods  to  validate  the  physical  parameters  of  a 
stripline  and  investigate  its  characteristic  impedance,  as  mentioned  previously.  This 
thesis  focuses  on  an  integral  equation  formulation  based  on  a  parallel  plate  Green’s 
function  and  efficient  entire  domain  method  of  moments  solution  for  the  center  strip 
current.  The  ground  planes  are  modeled  as  perfect  electric  conductors  (PECs).  Also, 
as  stated  in  Section  1.4,  the  center  conductor  is  modeled  as  infinitesimally  thin  and 
the  structure  is  considered  invariant  along  the  guiding  axis. 

1 . 6  Organization 

Chapter  2  provides  the  necessary  theoretic  background  starting  with  an  overview 
of  transmission  line  theory,  stripline  design  and  applications.  Chapter  3  provides  the 
development  of  the  required  parallel  plate  vector  potential  Green’s  function.  Chap¬ 
ter  4  details  the  formulation  of  a  generalized  pair  of  coupled  electric  field  integral 
equations  (EFIE’s).  It  includes  the  development  of  a  specialization  for  the  transverse 
electromagnetic  (TEM)  propagating  mode.  Chapter  5  details  the  Method  of  Moments 
(MoM)  solution  for  the  general  and  specialized  integral  equation  cases.  It  will  include 
a  discussion  on  the  type  of  expansion  and  testing  functions  selected  and  conclude 
with  numerical  results.  In  Chapter  6,  the  characteristic  impedance  of  the  stripline 
will  be  investigated  for  various  center  strip  alignment  configurations.  Proposed  design 
parameters  will  be  validated  using  this  information.  Lastly,  the  thesis  concludes  in 
Chapter  7  with  a  summary  of  research  efforts,  results,  potential  contributions  and 
future  research  opportunities. 
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II.  Background 


The  purpose  of  this  chapter  is  to  provide  a  background  for  the  stripline  field 
applicator,  and  to  explore  current  research  efforts  in  this  topic,  ft  is  assumed 
that  the  reader  is  familiar  with  electromagnetic  theory  and  the  uses  and  application 
of  Maxwell’s  equations.  This  chapter  will  look  at  basic  transmission  line  theory,  basic 
stripline  characteristics  and  design  practices,  applications  in  material  characterization 
and  current  efforts  to  mitigate  characterization  error  resulting  from  characteristic 
impedance  uncertainty. 

2.1  Transmission  Line  Theory 

The  stripline  belongs  to  a  family  of  devices  called  transmission  lines.  A  trans¬ 
mission  line  is  a  structure  that  forms  all  or  part  of  a  path  from  one  place  to  another 
for  directing  the  transmission  of  energy,  such  as  electromagnetic  waves  or  acoustic 
waves,  as  well  as  electric  power  transmission  [31].  Moreover,  in  electromagnetics,  a 
transmission  line  is  comprised  of  conductors  linking  electrical  systems  that  guide  en¬ 
ergy.  Additional  examples  of  transmission  lines  include  wires,  coaxial  cables,  electric 
power  lines,  and  microstrips.  Transmission  line  theory  must  be  used  when  the  length 
of  the  devices  to  be  analyzed  is  long  with  respect  to  the  operational  wavelength.  This 
relatively  large  size  allows  for  variation  in  voltage  and  current  distributed  throughout 
the  system  or  circuit. 

A  distributed  model  for  a  section  of  transmissions  line  is  shown  in  the  Figure 
2.1  where  R  is  the  resistance  in  both  conductors  per  unit  length  (f 1/m),  L  is  the 
inductance  in  both  conductors  per  unit  length  ( H/m ),  G  is  the  conductance  of  the 
dielectric  media  per  unit  length  (S/m)  and  C  is  the  capacitance  between  the  con¬ 
ductors  per  unit  length  ( F/m ).  [24]  R  is  included  to  account  for  ohmic  losses  in  the 
line.  G  is  included  to  account  for  dielectric  losses  in  the  dielectric  between  conduc¬ 
tors.  L  and  C  are  included  to  account  for  the  electric  and  magnetic  energy  stored 
in  the  line  [23].  For  lossless  lines  R  and  G  are  zero.  This  model  can  be  used  to  find 
fundamental  transmission  line  properties,  namely  the  propagation  constant  and  the 
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Figure  2.1:  Model  for  transmission  line  using  circuit  compo¬ 
nents. 


characteristic  impedance.  The  propagation  constant  ,  7  =  a  +  j/3,  is  calculated  in 
terms  of  the  inductance,  capacitance,  conductance  and  resistance  of  the  line  and  is 
given  by  the  expression  [23] 


If  the  line  is  lossless  (R,G  =  0),  the  propagating  constant  is  purely  imaginary.  For 
TEM  modes,  the  wave  number,  k,  is  related  to  the  propagating  constant  according 
to  the  following  relation 

7  =  jk  (2.2) 

where 

k  =  f3  —  ja  (2.3) 


and  behaves  like 


e-72  _  g -jkz  _  e-aze~jf3z 


(2.4) 
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with  a  =  0  for  a  lossless  line.  The  characteristic  impedance  of  a  transmission  lines  is 
the  ratio  of  the  voltage  to  the  current  for  a  single  traveling  wave.  The  characteristic 
impedance  has  units  of  ohms  and  is  given  as 


_  v_l_  _  [rJjZl 

°  If  \jc,+3<jC 


(2.5) 


When  using  transmission  lines  in  conjunction  with  other  devices  (i.e.  Genera¬ 
tors,  loads  or  other  transmission  lines)  the  behavior  of  the  energy  at  the  junctions  is 
important  to  understand.  For  maximum  power  transfer,  the  impedance  of  all  devices 
in  a  system  must  be  matched.  Variations  in  the  impedance  results  in  reflections.  The 
ratio  of  forward  traveling  energy  to  reflected  energy  is  known  as  the  reflection  coef¬ 
ficient.  The  reflection  coefficient,  expressed  in  terms  of  characteristic  impedance,  is 
given  by 


Z_l  ~  Z0 
ZL  +  Z o’ 


(2.6) 


where  Zl  is  the  load  impedance  and  Zq  is  the  characteristic  impedance  of  the  line. 


2. 2  Stripline 

Since  their  conception  in  the  1950’s  [31]  striplines  have  been  carefully  studied. 
Like  the  coax,  the  stripline  supports  a  dominant  transverse  electromagnetic  mode.  Its 
geometry  is  characterized  by  a  metal  strip  sandwiched  between  two  parallel  ground 
planes.  The  dominant  mode  of  the  stripline  is  non-dispersive  and  has  zero  cutoff 
frequency.  Striplines  are  also  relatively  easy  to  fabricate. 

Throughout  this  research,  the  following  coordinate  system  will  be  adopted  to 
describe  the  stripline  geometry.  A  cross  section  is  illustrated  in  the  Figure  2.2.  The 
origin  is  placed  at  the  center  between  the  ground  planes.  The  perfect  electric  con¬ 
ductor  (PEC)  ground  planes  located  at  y  =  ±h  extend  infinitely  in  the  x  and  z 
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Figure  2.2:  Stripline  cross  section  with  misaligned  center  con¬ 
ductor . 

directions.  The  infinitesimally  thin  center  conductor  is  located  at  y  =  ya  and  extends 
from  —a  <  x  <  a. 

2.2.1  Asymmetric  Striplines.  A  second  class  of  stripline  has  begun  to 
emerge.  This  class,  referred  to  as  asymmetric  striplines,  incorporates  all  stripline  that 
have  a  non-uniform  center  conductor  or  a  center  conductor  that  is  not  perfectly  cen¬ 
tered  between  the  ground  planes.  Asymmetry  effects  the  stripline  properties  presented 
in  previous  sections.  The  effects  of  the  misalignment  on  characteristic  impedance  will 
be  explored  in  Chapter  6. 

2. 3  Stripline  Design 

The  main  considerations  for  stripline  design  are  selecting  materials  used  in  fabri¬ 
cation,  calculating  physical  parameters  and  devising  a  method  to  feed  the  apparatus. 
Other  key  factors  are  operating  frequency  regime  and  characteristic  impedance.  In 
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Table  2.1:  Property  summary  of  material  for  center 
conductor  and  ground  plane  fabrication. 


Metal 

Conductivity 

Silver 

6.3  x  107 

Copper 

5.9  x  107 

Gold 

4.5  x  107 

Aluminum 

3.7  x  107 

Tin 

9.17  x  106 

H-214*(  at  2200  deg  F) 

8.1  x  105 

H-230*(at  1800  deg  F) 

7.9  x  107 

HR-120*(at  2200  deg  F) 

7.7  x  105 

H-282*(at  1800  deg  F) 

7.6  x  107 

*  High  temperature  alloys  produced  by  Haynes  International,  Inc. 

this  section  some  of  the  possible  fabrication  materials  will  be  reviewed  and  the  method 
of  calculating  the  physical  characteristics  of  the  stripline  will  be  described. 

Material  selection  must  begin  early  in  the  design  process  [20] .  Careful  selection 
of  the  dielectric  material  between  the  ground  planes  and  the  center  strip  will  strongly 
effect  the  performance  of  the  stripline.  Material  uniformity,  homogeneity,  isotropy, 
useful  temperature  range,  durability  and  frequency  response  are  the  main  factors 
governing  material  selection.  Table  2.1  lists  possible  choices  for  the  conductors. 

Physical  parameters  (such  as  operational  frequency  and  Z o)  are  governed  by  the 
intended  use  and  are  typically  user  defined.  The  first  step  in  specifying  the  operational 
frequency  regime  of  a  stripline  is  to  calculate  the  spacing  of  the  ground  planes.  The 
fact  that  the  ground  planes  form  a  parallel  plate  waveguide  will  be  exploited.  To 
ensure  that  no  higher  order  modes  are  excited,  the  cut-off  frequency  for  the  first 
higher  order  mode  above  the  dominant  mode  must  be  set  (as  a  rule  of  thumb)  ten 
percent  above  the  desired  highest  frequency  of  operation.  Next,  the  width  of  the 
center  conductor  must  be  computed  to  achieve  the  desired  characteristic  impedance. 
The  secondary  dimensions  of  the  ground  planes,  length  and  width  must  be  selected 
according  to  the  intended  application. 
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For  a  stripline  that  will  operate  in  the  dominant  mode  up  to  4  GHz,  stripline 
spacing  can  be  calculated  using  the  following  well  known  expression  [11] 

2fc  =  d-  (2-7) 

^ JCOl 

where  2 h  is  the  spacing  between  the  ground  planes,  c  is  the  speed  of  light  in  a  vacuum 
and  fCOl  is  the  cutoff  frequency  of  the  first  mode  above  TEM. 

An  approximate  value  for  the  center  conductor  width  can  be  found  by  using  the 
predetermined  center  conductor  width  (a)  to  the  ground  plane  spacing  ( h )  ratio  given 
by  the  expression  [10] 

1  =  ^  <2-8> 

More  exact  results  for  the  width  of  the  center  conductor  will  be  computed  using  a  full 
wave  analysis  in  Chapter  6. 

The  stripline  applicator  will  be  connected  to  a  network  analyzer  by  coax  cable 
during  operation.  Selection  of  a  feed  method  and  connector  type  can  effect  the  balance 
in  the  system.  Mismatches  between  the  characteristic  impedance  of  the  stripline  and 
coax  cable  will  introduce  error  into  material  characterization  measurements.  In  the 
Section  2.5  efforts  to  improve  junction  design  will  be  explored. 

2.4  Applications 

Striplines  have  many  applications.  They  have  becomes  staples  in  the  test  com¬ 
munity.  Striplincs  are  used  to  measure  the  electromagnetic  properties  of  materials 
and  to  test  the  radiated  immunity  of  automobile  components.  This  section  briefly 
summarizes  these  applications.  The  next  section  will  discuss  efforts  to  improve  the 
stripline  junction  design. 

2-4-1  Material  Measurements.  Striplines  are  frequently  employed  to  deter¬ 
mine  the  electromagnetic  properties  of  materials.  [9,39]  The  advent  of  modern  network 
analyzers  has  made  automated  calculation  of  the  constituent  parameters  possible. 
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Accurate  extraction  of  the  permittivity  and  permeability  of  a  material  require  a  well 
calibrated  system  and  reliable  extraction  methods. 


2.4. 1-1  Calibration.  One  method  of  calibrating  a  2-port  networks  is 
composed  of  through,  reflect  and  line  (TRL)  measurements  [9] .  The  thru  measurement 
is  accomplished  by  connecting  the  ends  of  the  port-1  and  port-2  terminals.  The  reflect 
measurement  is  accomplished  by  connecting  a  short,  or  reflecting  standard,  to  each 
port.  For  the  line  measurement,  a  length  of  transmission  line  is  inserted  between  the 
test  ports.  This  method  is  very  useful  in  waveguide  measurements. 

Other  test  beds  pose  a  challenge  to  calibration.  The  TRL  method  above  doesn’t 
lend  itself  well  to  the  stripline.  For  stripline  calibration,  an  empty  line  measurement 
and  three  carefully  placed  shorts  provide  the  information  necessary  for  calibration 
[9,39], 


2-4- 1-2  Extracting  Constitutive  Parameters.  The  Nicolson-Ross-Wier 
algorithm  (NRW)  is  used  to  calculate  the  permittivity  and  permeability  of  unknown 
materials  from  forward  scattering  parameters,  Sll  and  S21,  or  reverse  scattering  pa¬ 
rameters,  S22  and  S12  [9].  It  is  limited  to  single  layered,  homogeneous,  isotropic 
materials.  For  layered  material,  it  is  necessary  to  convert  the  scattering  parameters 
to  wave  transmission  parameters  (A-parameters)  and  isolate  the  A-parameters  of  the 
unknown  layer.  Once  isolated,  the  A-parameters  of  the  unknown  layer  can  be  con¬ 
verted  to  S-parameters  of  the  unknown  and  NRW  can  be  used  to  find  the  constituent 
parameters  for  the  unknown  layer  [9,39]. 

2-4-2  Automotive  Industries  Use  of  the  Tri-Plate.  The  automobile  industry 
uses  a  stripline-like  device,  a  tri-plate,  to  test  the  radiated  immunity  of  automobile 
components.  The  tri-plate  generates  an  electric  field  between  the  center  plate  and  the 
outer  plates.  As  automotive  components  are  place  in  the  tri-plate,  the  test  exposes 
both  the  component  under  test  and  its  wiring  harness  to  the  electric  held.  The 
operation  of  the  device  is  monitored  as  it  is  exposed  to  these  fields  [15]. 
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2.5  Efforts  to  Optimize  Striplines 

Mitigating  the  effects  of  junction  areas  has  been  a  major  thrust  in  research,  as 
indicated  by  the  many  reference  available  [6,15,16,27,30,32,34],  A  common  form  of 
coaxial-to-stripline  transition  consists  of  a  simple  inline  butt  joint  [26,32],  This  junc¬ 
tion  is  not  perfectly  matched  due  to  the  disparity  in  the  dimensions  of  the  two  pieces 
of  transmission  line.  Two  approaches  to  improve  impedance  matching  have  been  iden¬ 
tified  in  literature,  a  tapered  center  conductor  and  a  taper  in  conjunction  with  tuning 
wedges.  Similar  efforts  have  been  undertaken  by  the  Automobile  industries  to  perfect 
the  tri-plate  line  [15,27,30], 

In  recounting  the  design  of  his  measurement  device,  Barry  goes  into  great  detail 
about  his  efforts  to  design  the  ’’best  possible  transition”  [6],  The  optimal  junctions 
he  proposed  included  a  twenty  degree  taper.  These  tapers  were  introduced  to  avoid 
capacitances  between  the  conducting  strip  and  conducting  plane  where  the  connectors 
are  affixed  [34],  Barry’s  efforts  to  improve  the  junction  reduced  reflections  by  more 
than  30  dB  [6], 

Later,  Barry’s  method  was  also  combined  with  tuning  wedges  to  further  improve 
impedance  matching  by  Hanson  [16],  In  this  instance  an  eighteen  degree  taper  was 
used.  Tuning  wedges  were  made  by  inserting  tapered  corners  at  the  ends  of  the 
ground  plane  above  and  below  the  center  conductor.  However,  Hanson  concluded 
that  impedance  matching  was  not  critical  if  the  correct  method  was  used  to  calibrate 
the  system  and  extract  the  constitutive  parameters.  Non-linear  tapers  and  wedges 
composed  of  arbitrary  permittivity  profile  have  also  been  used  with  good  results  [37]. 

Tri-plates  incorporate  a  tapered  center  ground  plane  and  a  wedge  shaped  tran¬ 
sition  area.  Hwang  carefully  analyzed  the  impact  that  the  length  of  the  transition 
zone  and  the  steepness  of  the  wedge  have  on  energy  transmission  [21],  He  determined 
that  at  low  frequencies,  reflections  are  dominated  by  the  difference  in  height  between 
the  feed  and  test  region  of  the  apparatus.  He  proposed  that  a  gradual  discrete  taper 
would  provide  the  best  matching. 
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2. 6  Summary 

In  this  chapter,  fundament  concepts  from  transmission  theory  were  summarized. 
The  properties  of  striplines  were  examined.  Additionally,  the  stripline  design  process 
and  a  few  stripline  applications  were  presented.  The  chapter  concluded  with  a  brief 
summary  of  recent  developments  in  stripline  fabrication.  The  information  in  this 
chapter  will  be  used  throughout  the  rest  of  this  thesis.  The  geometry  described  in 
Figure  2.2  will  be  used  in  the  next  chapter  to  develop  the  parallel  plate  Green’s  func¬ 
tion.  The  approximate  values  for  the  ground  plane  spacing  and  the  center  conductor 
width  will  be  refined  in  Chapter  6. 
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III.  Parallel— Plate  Green’s  Function  Development 

3. 1  Introduction 

This  chapter  summarizes  the  derivation  of  the  vector  potential  Green’s  func¬ 
tion  for  an  electric  current  immersed  in  a  PEC  parallel-plate  waveguide,  as  seen  in 
Figure  3.1.  The  derivation  draws  heavily  upon  the  work  in  [17,22], 

The  vector  potential  parallel  plate  Green’s  function  is  derived  by  first  finding 
the  boundary  conditions  in  terms  of  vector  potential.  Next,  the  general  solution  of  the 
vector  potential  wave  equation  is  found.  Finally,  boundary  conditions  are  enforced  on 
the  general  solution  leading  to  a  unique  representation  of  the  vector  potential  in  terms 
of  the  Green’s  function  and  the  electric  current.  Each  of  these  steps  is  summarized 
below. 

3.2  Vector  Potential  Boundary  Conditions 

Enforcing  boundary  conditions  (BCs)  on  the  electric  field  (3.1)  is  complicated 
by  the  required  gradient  and  divergence  operators  [5], 

E  =  — - — (k2  +  VV-)A  (3.1) 

ju/de 

However,  these  operations  can  be  avoided  by  recasting  the  BCs  in  terms  of  the  vector 
potential.  This  is  done  by  taking  advantage  of  the  non-coupling  nature  of  a  current 
in  a  planar  PEC  environment  (  i.e.,  each  component  of  current  maintains  only  its 
respective  component  of  potential)  [17].  The  vector  potential  BCs  are  found  to  be 


Ax(x,y  =  ±h,z)  =  0 

V  x,z 

(3.2) 

dAy(x,y  =  ±h,z) 
dy 

V  x,z 

(3.3) 

Az(x,  y  =  ±h,z)  =  0 

V  x,z 

(3.4) 

Since  the  BC’s  presented  in  equations  (3.2)-(3.4)  must  hold  for  all  values  of  x, 
z  (i.e.,  —  oo  <  x  <  oo,  —  oo  <  z  <  oo),  a  Fourier  transformation  on  these  spatial 
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Figure  3.1:  Problem  Geometry:  General  3D  Current  (J)  Im¬ 
mersed  in  Parallel  Plate  Waveguide. 

variables  is  performed.  Transforms  will  be  performed  according  to  the  defined  forms 
below  (see  equations  (3.5)  and  (3.6)  ).  Spectral  functions  are  indicated  by  a  tilde 
overbar  and  functional  dependance  on  spectral  variables  (x  — »  £,  y  — >  rj  and  z  — ■»  (). 
A  double  transform  is  equivalent  to  two  single  transforms. 


/OO  1  POO 

f(x,  y,  z)e~Jz(dz  f(x,  y,z)  =  —  /  f(x,  y,  ()eJzCd(  (3.5) 
-OO  J  —  OO 

~  POO  1  POO  _ 

f(Z,y,  0  =  /  f(x,y,()e~33*dx&  f(x,y,  ()  =  —  /(£,  2/,  (3.6) 

Transforming  the  vector  potential  boundary  conditions  to  the  £,  domain  leads 
to  the  following  spectral  representation  of  the  vector  potential  boundary  conditions, 


Ax{£,y  =  ±h,  C)  =  o 


(3.7) 
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dAy(^,y  =  ±h,  C)  _ 

dy 

Az(£,y  =  ±h,C )  =  0 

As  in  equations  (3.2)-(3.4),  equations  (3.7)-(3.9)  indicate  that  the  tangential  compo¬ 
nents  of  the  vector  potential  vanishes  at  the  PEC  boundary;  similarly,  the  normal 
derivative  of  the  normal  component  is  also  zero.  These  boundary  conditions  will  be 
used  to  find  the  desired  Green’s  function.  The  general  solution  to  the  vector  potential 
wave  equation  is  found  next. 

3.3  Solution  of  the  Vector  Potential  Wave  Equation 

A  general  solution  to  the  vector  potential  wave  equation, 

V2A(x,  y,  z)  +  k2A(x ,  y,  z))  =  -yJ( \x,  y,  z),  (3.10) 

can  be  found  as  the  linear  superposition  of  principal  and  reflected  contributions.  The 
principal  contribution  is  related  with  a  general  current  immersed  in  a  homogenous 
space  similar  the  the  material  found  between  the  ground  planes  of  the  parallel  plate 
waveguide.  The  reflected  contribution  is  due  to  the  presence  of  the  ground  planes. 
The  respective  contributions  satisfy  the  differential  equations  given  by 

V2Apa(x,  y,  z)  +  k2Apa(x ,  y,  z)  =  -yJa(x,  y,  z) . . .  a  =  x,y,  z  (3.11) 

V2Ara(x,  y,  z)  +  k2Ara(x ,  y,z)  =  0  . .  .a  =  x,y,z  (3.12) 

where  the  superscript  p  denotes  the  principal  contribution  of  the  general  wave  equation 
and  the  superscript  r  denotes  the  reflected  contribution  of  the  general  wave  equation. 
The  total  solution  is 


(3.8) 

(3.9) 


4>(z,  V,  *)  =  A*(x,  y,  z)  +  Ara(x,y,z)...a  =  x,  y,  z  (3.13) 
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The  solutions  to  the  principle  and  the  reflected  parts  are  described  in  the  sections 
that  follow. 


3.3.1  Principal  Solution.  The  principal  part  of  the  vector  potential  wave 
equation(3.11)  can  be  solved  using  a  Fourier  transform  based  complex-plane  analysis. 
Transforming  on  z  and  x  according  to  the  transform  pairs  (3.5)  and  (3.6) (with  the 
aid  of  the  Fourier  differentiation  theorem)  leads  to 


d2APa(£,y,  C)  2  7P  (  \  1 

- Q~2 - P  Aa(x ,  Vi  z )  =  ,  y,0---ot  =  x,y,z  (3.14) 

where  p  =  a/£2  +  C2  —  k2.  Note,  £  is  an  x-directed  wave  number  playing  the  role  of 
kx,  p  is  a  ^/-directed  wave  number  playing  the  role  of  jky  and  (  is  a  ^-directed  wave 
number  playing  the  role  of  kz. 

Since  ff^.(x,r/,  z)  is  associated  with  a  current  in  an  unbounded  region,  a  Fourier 
transform  on  y  can  also  be  performed  according  to  the  definition  below,  namely 


0 


f~(Z,y,Oe-jvydy^ht,y,() 


^  /  m,yX)ejvydy. 


(3.15) 


Completing  the  desired  transform  leads  to 


-{rj2  +p2)Apa(£,ri,() 
k(Z,V,  0 


PJQ(£,V,  0 
(rj2  +  p2) 

pjg(t,y,  0 

(v  +  jp)(v-jp) 


(3.16) 


Since  boundary  conditions  will  need  to  be  enforced  on  the  total  solution  at 
y  =  ±h,  an  inverse  Fourier  transform  must  be  performed  on  the  principal  contribution 
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Figure  3.2:  Singularities  of  the  Integrand  of  Equation  (3.17) 
in  Complex  rj  Plane. 

to  return  it  to  the  ^/-domain.  The  inverse  transform, 

OO  OO  ~ 

*  0  =  kb*  <W*n  =  ^  (3'17> 

— OO  — OO 

is  accomplished  with  the  aid  of  complex  analysis,  as  discussed  next. 

3. 3. 1.1  Singularity  Location  and  Closure  Conditions.  For  the  inte¬ 
grand  of  equation  (3.17),  the  only  singularities  are  simple  poles,  rj  =  ±jp.  The  loca¬ 
tion  of  these  poles  can  be  determined  by  examining  the  expression  p  =  \/£2  +  C2  —  k2. 
Note,  £  and  £  are  the  integration  variables  along  the  real  axis  and  are  purely  real. 
Furthermore,  for  passive  media  k,  given  by  k  —  will  have  a  positive  real  part 

and  a  negative  imaginary  part.  Hence,  k 2  will  be  found  in  the  third  and  fourth  quad- 
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rants  of  the  complex  plane.  Likewise,  p 2  ( p 2  ~  —A;2)  will  be  found  in  the  first  and 
second  quadrants  of  the  complex  plane.  Taking  the  root  of  p2,  p  and  —  p  will  be 
found  in  the  first  and  third  quadrants  of  the  complex  plane  respectively.  Finally,  the 
multiplication  of  the  value  of  p  with  j  will  rotate  the  location  of  the  pole  into  the 
second  and  forth  quadrant  of  the  complex  //-plane  (See  Figure  3.2). 


Cauchy’s  integral  theorem  can  be  used  in  the  evaluation  on  the  real-axis  integral 
in  equation  (3.17).  An  appropriate  contour  for  the  case  of  y  >  y'  is  shown  in  Figure 
3.3.  Note,  since 


ejv(y-y')  —  e-Vim(y-y')  ejvre(y-y’) 


(3.18) 


a  careful  inspection  reveals  that  the  following  closure  conditions  must  be  obeyed  for 
convergence  [17]: 


1.  Case  y  —  y'  >  0:  e~‘nim^y~y'^  provides  decay  if  yim  >  0  and  an  upper-half  plane 
closure  (UHPC)  will  ensure  the  infinite  semi-circular  contribution  vanishes  as  y 
approaches  oo. 

2.  Case  y  —  y'  <  0:  e~mrn<<y~^  provides  decay  if  yim  <  0  and  a  lower-half  plane 
closure  (LHPC)  is  required  to  ensure  the  infinite  semi-circular  contribution  van¬ 
ishes  as  y  approaches  oo. 

Thus,  the  relation  of  the  observer  location  ( y )  to  the  source  location  (//')  deter¬ 
mines  whether  an  LIHPC  ( y  >  y')  or  a  LHPC  ( y  <  y')  is  required. 

3.3. 1.2  Applying  Cauchy  Integral  Formula.  Now  that  the  singularities 
have  been  identified  and  the  closure  conditions  investigated,  we  may  now  proceed  in 
finding  A^,(£,  y,  ().  Only  the  case  of  y  >  y'  will  be  presented  here. 

Recall  that  Ja(£,  y,  ()  can  be  written  according  to  the  following  Fourier  trans¬ 
form  pair, 
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Figure  3.3:  Upper-Half  Plane  Closure  used  to  integrate  case 

y>y' 


~  poo  „  ph  _ 

Mt,V,C)=  UWX)e-^'dy'=  ja{U,Qe-jvy'dy'.  (3.19) 

J—  oo  J  —  h 

Note,  the  region  of  existence  of  the  current  is  physically  confined  within  the  parallel- 
plate  background  environment  — h  <  y  <  h ;  accordingly,  the  limits  of  integration  can 
be  truncated  to  neglect  the  portions  beyond  the  boundaries  that  don’t  contribute  to 
the  integral  (i.e. ,  Ja(x,h  <  y  <  oo,  z)  =  0)).  Substituting  (3.19)  into  (3.17)  leads  to 


K&v,  0  =  ^ 


M  f  M£,y',()e  jny' dy' 

—h 


(v  +  jp)(v~  jp) 


-ejvydr] 


eiv{y-y') 


2vr  ,/  {y  +  jp)(y-jp) 

—h  \  — oo 


—dr]  fj,ja(£,y',  Qdy'  (3.20) 


where 


&a(tX-,y-y') 


l  r  e  jv(y-y') 

2t t]  (  f]  +  jp)  (rj  -  jp) dn 

—  OO 


(3.21) 
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is  the  spectral  domain  (i.e.,  the  £,  (-  domain)  Green’s  function  for  the  principal  wave 
contribution. 

A  closed  form  expression  for  the  spectral  domains  Green’s  function  is  calculated 
next  using  complex  plane  analysis.  Only  the  process  for  the  case  y  >  y'  (UPHC)  will 
be  presented  (see  Figure  3.3).  Momentarily  neglecting  the  factor  A  in  the  spectral  do¬ 
main  Green’s  function,  the  desired  inverse  Fourier  transform  integral  to  be  calculated 

is 


em{y-y ') 


drj  =  lim 


ejv(y-y') 


(y+jp)(y~jp)  R^oo  J  (rj  +  jp)(r)  -  jp) 

Cre 


—,d'n 


R 


=  lim 


ojy{y-y') 


R^ooJ  (r]  +  jp)(r)  -  jp) 

-R 


drj 


(3.22) 


Applying  Cauchy’s  Integral  Theorem  [3]  leads  to 


R 


ojy(y-y') 


ifco  J  (rj  +  jp)  (rj  -  jp) 


rrd'//  + 


Jy(y-y') 


-R 


c+ 


(v  +  jp)(v-jp) 


drj 


ejy(y-y') 


lim  /  - - — - -dp  =  0 

r^oc  J  (rj  +  jp)(rj  -  jp) 

c+ 


R 


Jv(y-y') 


J  (rj  +  jp)  (rj  -  jp) 

—R 

e  jy{y-y') 


•  \drj  = 


cjy(y-y') 


(3.23) 


— d/rj  —  lim 


c+ 


(y  +  jp)  (rj-  jp)  R->oo  J  (rj  +  jp)(r)  -  jp) 

c+ 


—dp  (3.24) 


The  left-hand  side  of  equation  (3.24)  represents  the  desired  integration  along  the  real 
axis.  The  right-hand  side  of  equation  (3.24)  is  solved  using  Cauchy’s  integral  formula 
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and  Jordan’s  Lemma  [3],  leading  to 


~§ 


ojvi.y-y') 


c+  (v  +  3P)(v~3P) 

ejv{y-y') 


dr]  =  — 


-j  2vr 


(V  +  3P) 


(~i-X 


\n =jp 


eMy-v') +  jp) 

(v  ~  JP) 

_  ne-p{y-y,) 

P 


dr] 


(3.25) 


It  should  be  noted  that  the  negative  of  the  residue  is  used  to  compensate  for  the 
clockwise  direction  of  integration.  The  infinite  semicircular  contribution  (C+)  van¬ 
ishes  according  to  Jordan’s  lemma.  Furthermore,  application  of  Jordan’s  lemma  can 
be  extended  to  y  —  y',  as  the  integrand  still  converges  as  strongly  as  \  for  large  y. 

The  combined  result  of  both  cases  y  >  y'  and  y  <  y'  is  (the  factor  of  \/2n  again 
once  again  included) 

~  P-p\y-y'\ 

GL.K.C;  !/-!/)  =  — 5— -  (3.26) 


which  leads  to 


2  p 


/e-p\y-y'\  ~  r  e-v\y-y'\  ~ 

— 2^ — pJady'  =  j  — — — ]iJady 
y*  -h 


(3.27) 


3.3.2  Reflected  Solution.  The  form  of  the  reflected  portion  of  the  vector 
wave  equation  in  the  £,  ("-domain  is  easily  determined 


y,Q 

dy2 


~  P  Ara(£,  y,  C)  =  0  . . .  a  =  x,  y,  z. 


(3.28) 


Its  harmonic  solution  is  readily  recognized  to  be,  Ara  =  W^e~py  +  W~epy.  Note,  the 
first  term  represents  an  up-going  wave  with  complex  amplitude  Wfl  and  the  second 
term  represents  a  down-going  wave  with  complex  amplitude  W~. 


3-4  Applying  Boundary  conditions  to  Total  Solution 

Applying  the  boundary  conditions;  (3.7)-(3.9)  to  the  total  solution,  Aa  =  Ap  + 
Ara ,  at  y  =  ±h  yields 
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An  — 


o-p[y-y'\ 


2  p 


-pJdy'  + 


f>  p-p(y'+2h+y)  ,  z~)  D  p-p{-y'  +lh+y) 

1Laac'  '  JLaa-LLaac' 


2p(l  -  RaaRaae~4Ph) 


f>  p-p(~y'+2h-y )  id  p  p—p(y'+ih—y 
iiqqg  i  itaaJtaQ!c' 

2p(l  -  RaaRaae~Aph) 


pJady' 


(3.29) 


where  Raa  is  the  reflection  coefficient  for  the  lower  ground  plane,  Raa  is  the  reflection 
coefficient  for  the  upper  ground  plane  and  p  =  \/£2  +  (2  —  k2.  All  of  the  integrand  in 
equation(3.29)  being  multiplied  by  pJ  represents  the  vector  potential  parallel-plate 
Green’s  function  in  the  spectral  (£,  ()  domain. 

In  this  form  (3.29),  a  great  deal  of  physical  meaning  can  be  gleaned.  As  seen  in 
Figure  3.4,  each  term  in  the  solution  accounts  for  a  different  possible  path  from  the 
source,  y'  ,to  the  observer,  y  [17]. 


•  The  1st  path,  direct  from  source  to  observe,  is  accounted  for  in  the  contribution 
from  the  principal  part,  e~p^v~y''1 

•  The  2nd  path,  from  the  source  to  the  observe  after  reflection  off  the  upper- 
conductor,  is  accounted  for  in  the  term  Re~p^~y' +2h~y'> 

•  The  3rd  path,  from  the  source  to  the  observer  after  reflection  off  the  lower- 
conductor,  is  accounted  for  in  the  term  Re~ptyV' +2/HA) 

•  The  4th  path,  from  the  source  off  the  top-conductor,  then  the  bottom-conductor 
then  to  the  source,  is  accounted  for  in  the  term  RRe~p^~v> +ih+y) 

•  The  final  path,  the  inverse  of  the  4th  path,  is  accounted  for  in  the  5th  term 
RRe-p(y'+4h-y) 

Clearly,  the  proposed  principal  contribution  is  related  to  the  forced  response  in 
unbounded  media  while  the  reflected  contribution  accounts  for  the  reflected  response. 
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o  e  ©  ©  o 
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y  — 

v  — 


-h  — 

Figure  3.4:  Each  Possible  path  (1-5)  from  the  source  to  the 
observer  is  accounted  for  in  equation  (3.29). 
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A  more  compact  version  of  the  Green’s  function  can  be  found  by  using  the 
definition  of  hyperbolic  functions  and  setting  Raa,  Raa  =  —  1 . . .  a  =  x,  z,  Ryy,  Ryy  = 
1.  Equation  (3.29)  reduces  to 


Gota(jL)  C;  y  I  Do) 


cosh\p(2h-\y-y'\)\-cosh\p(y+y')\  -r 
2psinh(2ph) 

cosh[p(2h-\y-y'\)]+cosh[p(y+y')]  -f  _ 

2psinh(2ph)  ” 


(3.30) 


Note,  the  Green’s  function  is  even  in  p,  thus  the  square  root  branch  cut  contribution 
is  removable,  as  physically  expected  for  this  closed  structure. 


3. 5  Summary 

In  this  chapter,  the  vector  potential  parallel-plate  Green’s  function  for  an  elec¬ 
tric  current  was  developed.  First,  spectral  and  spacial  vector  potential  boundary 
conditions  for  the  parallel-plate  wave  guide  were  found.  The  forced  vector  potential 
wave  equation  was  solved  by  decomposing  the  equation  into  principle  and  reflective 
parts.  Each  part  was  then  solved  independently.  Complex  analysis  was  used  to  find 
the  solution  to  the  principle  part  and  the  reflected  part  was  readily  identified  as  having 
a  harmonic  solution.  Boundary  conditions  were  applied  to  the  total  solution  at  the 
outer  conductors  to  find  the  Green’s  function.  The  hyperbolic  version  of  the  Green’s 
function  (3.30)  was  identified  and  will  be  used  in  the  next  chapter  to  formulate  an 
integral  equation  for  a  stripline. 
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IV.  Electric  Field  Integral  Equations 

4-1  Introduction 

In  chapter  3,  we  discussed  the  development  of  the  parallel-plate  Green’s  func¬ 
tion.  In  this  chapter,  we  will  use  this  Green’s  function  to  formulate  coupled  integral 
equations  for  hireling  the  current  on  the  center  conductor  of  a  stripline  modeled  as  a 
strip  embedded  in  a  parallel-plate  waveguide  environment.  First,  a  set  of  coupled  in¬ 
tegral  equations  will  be  formulated  for  a  strip  with  infinite  conductivity.  This  method 
will  include  converting  the  PEC  boundary  condition  to  the  ((-domain  (prompted  by 
the  infinite  length  along  the  z-axis).  Finally,  a  specialized  integral  equation  for  a 
stripline  supporting  a  TEM  mode  will  be  considered.  In  the  next  chapter,  the  method 
of  moments  will  be  used  to  solve  for  the  unknown  currents  and  guided  propagation 
constant  found  in  these  integral  equations. 

4-2  Integral  Equation  Formulation  for  PEC  Stripline 

The  goal  of  this  section  is  to  formulate  an  integral  equation  for  the  unknown 
current  density  on  the  PEC  center  conductor  of  the  stripline.  Specifically,  an  electric 
held  integral  equation  (EFIE)  is  developed  by  enforcing  the  well-known  boundary 
condition 

t  ■  E  —  0  . .  .t  —  x,  z  (4.1) 

on  the  PEC  center  strip.  Note,  it  will  be  assumed  that  the  PEC  center  strip  is 
infinitesimally  thin  and  located  at  y  =  y$. 

Decomposing  the  electric  field  in  (4.1)  into  incident  and  scattered  contributions, 
E  =  El  +  Es ,  leads  to  the  result  (for  t  —  x,z) 

t  ■  Es  =  -t  ■  El  or  Esx  =  -El  Esz  =  —Elz  (4.2) 

where 

Es  =  — (k2  +  VV-US.  (4.3) 

jujeg 
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Substituting  equation  (4.3)  into  the  leftside  of  the  equation  (4.2)  and  transform¬ 
ing  to  the  ((-domain  yields 


t  •  — (/c2  +  VVdA5  =  -t-  E{  ...t  =  x,z 

juefi 


(4.4) 


where 


V 

As 


,  d  „ d 
x—  +  y—+j(z 
ox  ay 


G( C;  p\p')'J (/£>',  Qdx'dy'  ...p  =  xx  +  yy...p 


y  x 


x'x  +  y'y ,  (4.5) 


G(p  |  p',(()  is  the  ((-domain  dyadic  Green’s  function,  y*,x*  is  shorthand  notation 
for  the  region  of  integration  and  J,  the  ((-domain  representation  of  a  the  £th  natural 
mode  z-propagating  surface  current  on  an  infinitesimally  thin  center  conductor,  can 
be  written  as  [28] 


f_  -  Jsx(x)S(y  -  y0)  ,  ,  Jsz(x)S(y  -  y0) 

'  J  —  00  /  ^  ^  N  |  Z 


(C  -  0) 


(C  -  0) 


(4.6) 


Additionally,  the  limits  of  integration  for  equation  (4.5)  can  be  reduced  to  the  region 
where  the  strip  exists  (—a  <  x'  <  a). 
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Substituting  equations  (4.5)  and  (4.6)  into  equation  (4.4)  yields  the  following 
boundary  condition  relation  for  the  x  and  z  components  of  the  electric  field 

1  _  f)2  ~  8 

- - (k2  Ax  +  +  —j(Az)  =  —  El 

juje/i  8xz  ox 


jLjefi 


[ik 2 
~2w 


&xx(S,C Wo  I  y0)ej^e-^^^-didx' 


—a  —  oo 
a  oo 


-  |y  /  /  ^Gxx^X'iVo  I  y, 


,)ejx^e  ix'iJ*Ar  ]  d^d.,-' 

(C  -  Ce) 


—a  — oo 
a  oo 


—a  —  oo 


~J  J  C  ;yo\Vo)e^e-^^^-)didx' 


(k2  -  £2)(L(£,C Wo  I  ?/0)e^e-^,?Js;c(x/)c/^/ 
C £&zz(Z,C-,yo  | 


=  -e. 


=  —2njue((  -  Ct)K 


(4.7) 


—  (k2Az  +  j(Ax  -  C2AZ ))  = 

jcueyu  ox 


juje/i 


fik2 

~2 V 


&UiX;  yo  I  Voy^e-^^AA-dtdx' 


—a  —  oo 
a  oo 


k  I  I  t rrt  it-  i  ..  \„jx£„—j x£  Jsx\x  ) 


2tt  y  y  1  y°)ejx  e  (c  -  Q) 


d^dx' 


—a  — oo 
a  oo 


~rj  j  |  *)«*«-*£!§ 


d£dx 


=  -El 


(£,Gxx(t,,C,yo  |  Vo)e3xie  Jx^J8X(x')d^dx' 

(k2  -  (2)dzz(^C -y0  I  y0)e3X^e~3X^  Jsz(x')d^ dx1 


—a  — oo 


=  -2ir^e(C  -  0)£” 


(4.8) 


For  natural  mode  propagation  (£  =  (^),  of  which  the  TEM  is  dominant,  the 
right-hand  side  of  equations  (4.7)  and  (4.8)  vanish  [28]  leading  to  the  desired  coupled 
EFIE’s  for  the  unknown  currents  Jsx  and  Jsz,  that  is 
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a  oo 


J  J(k 2  -  £2)<L(£,  C;  Vo  |  y0)e^e-^Jsx{x')didx' 

—  a  —  oo 

a  oo 

-  J  J  C&zz(S,  C;  Vo  |  y0)e^e-^  jsz(x')d&x'  =  0  (4.9) 

—a  — oo 

a  oo 

-  /  / C£<L(£,C;3/o  | 

—  a  — oo 

a  oo 

+  /  J  (k2  —  C2)4,(e,  C;  Vo  |  Vo)e?*e~i*  jsz^)d^dx'  =  0  (4.10) 

—a  — oo 


It  should  be  noted  that  the  first  integral  in  equation  (4.9)  accounts  for  the  x 
component  of  the  scattered  electric  held  maintained  by  the  x  component  of  current. 
The  second  integral  in  of  equation  (4.9)  accounts  for  the  x  component  of  the  scattered 
electric  held  maintained  by  the  z  component  of  current.  The  hrst  integral  in  equation 
(4.10)  accounts  for  the  z  component  of  the  scattered  electric  held  maintained  by  the 
x  component  of  current.  Finally,  the  last  line  in  integral  (4.10)  accounts  for  the  z 
component  of  the  scattered  electric  held  maintained  by  the  z  component  of  the  current. 
These  relationships  will  be  exploited  in  the  next  section  to  find  a  specialized  TEM 
integral  equation  and  will  subsequently  lead  to  enhanced  computational  efficiency. 

4-3  TEM  Specialization  of  a  PEC  Stripline 

For  TEM  mode  propagation  in  the  air  filled,  lossless  stripline,  there  is  no  x 
component  of  the  current  ,  that  is,  Jsx  =  0  and  the  propagation  constant  takes  on 
the  value  (  =  k  =  ko  (see  Section  5.3  for  more  information  on  these  results).  Under 
these  circumstances,  equation  (4.10)  vanishes  and  equation  (4.9)  reduces  to 

a  oo 

J  J  kot&zz(Z,  k0]  y0  |  yo)ej^e-^Jsz(x,)dfdx'  =  0  (4.11) 
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Equation  (4.11)  represents  a  single  EFIE  for  the  unknown  surface  current  Jsz  of 
the  Dominant  TEM  mode,  and  will  be  explored  in  the  next  chapter.  It  is  anticipated 
that  substantial  computational  gains  will  be  achieved,  compared  to  the  coupled  EFIE 
equations  (4.9)  and  (4.10).  Note,  if  the  properties  of  higher  order  modes  are  sought, 
the  coupled  EFIE  must  be  used. 

4  •  4  Summary 

In  this  chapter,  coupled  integral  equations  were  formulated  for  a  stripline  having 
a  perfect  center  conductor.  Furthermore,  it  was  shown  that  the  coupled  integral 
equation  reduce  to  a  single  integral  equation  for  the  TEM  mode.  This  specialization 
naturally  leads  to  a  more  computationally  efficient  formulation.  The  solution  of  these 
integral  equations  will  be  sought  in  the  following  chapter. 
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V.  Method  of  Moments 


5. 1  Introduction 

In  this  chapter,  an  efficient  entire-domain  method  of  moments  (MoM)  technique 
will  be  used  to  solve  the  integral  equations  formulated  in  the  last  chapter.  For  the 
coupled  integral  equations  formulated  for  the  perfectly  conducting  center  conductor, 
a  Galerkin’s  based  MoM  is  found  with  the  aid  of  Chebyshev  polynomials.  Solutions 
will  be  found  for  both  the  general  coupled  EFIE  formulation  and  the  specialized  TEM 
EFIE  formulation.  The  currents  ( Jsx  and  Jsz)  and  guided  propagation  constant  (£) 
found  for  each  of  these  formulations  through  this  process  will  be  used  in  the  next 
chapter  to  calculate  the  desired  characteristic  impedance  of  the  stripline. 

5.2  MoM  Solution  of  Coupled  EFIE’s  for  PEC  Strip 

In  this  section,  a  Galerkin’s  based  MoM  solution  is  sought  for  the  coupled  EFIE’s 
for  the  stripline  having  a  PEC  center  strip.  Two  different  scenarios  of  test  and  ex¬ 
pansion  functions  will  be  investigated.  It  will  be  shown  that  the  second  set  of  expan¬ 
sion  and  testing  functions  has  more  desirable  rates  of  convergence,  and  will  therefore 
be  used  in  solving  for  the  unknown  currents  and  calculation  of  the  characteristic 
impedance  and  TEM  mode  propagation  constant. 

5.2.1  General  MoM  Solution.  In  the  previous  chapter,  it  was  shown  that 
the  coupled  integral  equation  (CIE)  for  a  stripline  having  a  PEC  center  conductor  are 
given  by 
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a  /»oo 


—a  J  — oo 


a  roc 


a  /‘oo 


(k2  -  £2)(L(£,  C;  yo\yo)ejx^-^Jsx{x')didx' 

^zz{^yo\yo)ej^e-^Jsz{x')didx'  =  0 
£CGW£,  C;  yo\yo)ej^e-^Jsx{x')didx' 


—a  J  — oo 


—a  J  — oo 


a  /‘oo 


+ 


(&2  -  C2)4*(£,  C;  2/0|2/0)^e-^  Jaz{x')d£dx'  =  0. 


—a  J  —  oo 


(5.1) 

(5.2) 


The  first  step  in  the  MoM  solution  is  to  expand  the  unknown  currents  Jsx,  Jsz  accord¬ 
ing  to  the  weighted  sum  given  by 

N- 1 

Jsa(x')  =  ^2  anaena(x')  ...a  =  X,Z  (5.3) 

n= 0 

where  ena(x')  are  known  expansion  functions  and  ana  are  unknown  expansion  coeffi¬ 
cients.  Substituting  (5.3)  into  (5.1)  and  (5.2)  leads  to 


N—l 

n=0 


N-l 

E 

n=0 


(k2  -  £2)Gxx(£,C,yo\yo)e:>x!i  /  enx(x')e  Jx^dx'd^ 

)  J  —a 

/oo  _  pa 

^CG»K.C;S/o|l/o)e’*f  /  enI(x')e~”'(dx'd^  =  0 


(5.4) 


iV-1 


XX*  /  fCG**(f,C;yo|3/o)e,a*  /  ena:(a;')e  3X'^dx'd^ 


/  Qnx 
n= 0  ^  -°° 

N-l 


■' *  t-  oo  _  />a 

+  y,  an2  /  (/c2  -  C2)GZ*(£,  C;y0|?/0)eJX?  /  enz(x')e~JX'tdx'd£,  =  0.  (5.5) 


Q"riz 

n= 0  •/-°° 


Defining  .gnQ(£)  as 


,,(£)  =  /  enQ(a/)e  3X  ^dx' . . .  a  =  x,  z 


(5.6) 
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allows  (5.4)  and  (5.5)  to  be  simplified  as  follows 


N  1  /*oo  _ 

anx  /  ( k 2  -  £2)Gxx(£,  C;  yo\yo)e,x*gnx(€)dt 

n= 0  ^ 

N—  1  /»oo  ~ 

-  ^  Unz  /  =  0  (5‘7) 

n= 0  ^°° 

1  /*oo  ~ 

~^2anx  /  £CG*x(£, C; yo|yo)e^„x(0^ 

n=0  '/~°° 

1  /»oo  _ 

+  /  (fc2  -  C2)G^(e,C;?/o|?/o)e^(0^  =  0.  (5.8) 

Equations  (5.7)  and  (5.8)  represent  2  equations  and  2N  unknowns.  In  order  to  estab¬ 
lish  a  2N  x  2N  MoM  matrix,  the  following  testing  operators  are  applied  as  the  next 
step  in  the  MoM  solution,  namely 

a 

I1  IW*)*- ••“  =  *.*  and  m  —  0, . . .  ,N  —  1  (5.9) 

—a 

where  tma(x)  is  a  known  testing  function.  Applying  the  a  =  x  testing  operator  to 
(5.7),  the  a  =  z  testing  operator  to  (5.8)  and  defining  fma(£)  as 

a 

fma(€)  =  I  tma(x)eJX^dx  (5.10) 
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results  in  the  following  coupled  equation 


N—  1  poo 


£< 


(k‘  -  r)G„K,  C;  y„|y„)/™»K)9ra(0<*« 


^nx 

n= 0  ^  -°° 

N—l  poo  „ 

-  /  £C^zz(£>  C; yo|yo)/ms(0^n*(0d^  =  o . . .  m  =  o, . . . ,  n  -  1 

n=0  ^°° 


(5.11) 


1  /•oo  _ 

^  I  CCGxx(Ci  Ci  VolVo)  fmz(C)9nx(C)dC 

n=  0 

JV-1 


E' 


^nz  I 

n= 0  ^-°° 


(k  ~  C  )GZZ(C,  C;  Volvo)  fmz(C)gnz{C)d£  =  0...m  =  0,  ...,iV-l 


which  can  be  further  simplified  as 


where 


N 


N 


E'*--4”"  +  E“-4S"  =  0'"m=0'---"Ar-1 

n=0  n= 0 

AT  AT 

+  ^«n2bl™l  =  0...m  =  0,...,iV-l 


n=0 


n=0 


(5.12) 


(5.13) 


/j  mn  _ 


j^mn  _ 


j^rnn  _ 


a  mn  _ 


(k2  -  C2)GXX(C,  C;  yo\yo)fmx{C)9nx{C)dC 

OO  _ 

C(Gzz(C,  C;  Volvo)  fmx{C)9nz{C)dC 
C(Gxx(C,  C;  Volvo)  fmz(C)gnz(C)dC 


'  — oo 

POO 


J  — oo 
r»oo 


(■ k  -  C  )GZZ(C,C',  Volvo)  fmz(C)9nz(C)dC- 


(5.14) 

(5.15) 

(5.16) 

(5.17) 


Note,  A™  (with  a  and  (3  equal  to  x,  z )  physically  represents  how  the  nth  /3-directed 
source  current  term  couples  into  the  mth  CK-directed  field  term.  Next,  two  specific 
scenarios  of  expansion  and  testing  functions  will  be  investigated. 
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5.2.2  Galerkin  Method  Using  Chebsyshev  Polynomials  of  the  First  Kind. 
Selecting  the  ’’best”  expansion/testing  functions  has  been  shown  to  improve  the  ac¬ 
curacy  and  computational  efficiency  of  the  MoM  [4,5,19,29].  Using  an  entire-domain 
technique  can  further  improve  efficiency  because  fewer  terms  are  required  in  the  ex¬ 
pansion  leading  to  a  smaller  MoM  matrix  [5].  Chebsyshev  polynomials  make  excel¬ 
lent  expansion  and  testing  functions  in  this  research  because  of  their  ability  to  model 
current  behavior  on  a  strip,  include  edge  effects  [4, 14],  The  proposed  testing  and 
expansion  functions  in  the  first  scenario  are  given  by  [17] 


X  X 

tnx  =  enx  =  T,2n+ 1)(— )\/ 1  —  (-)2...(odd  function)  (5.18) 

a  V  a 


tnz  Cnz 


Wf) 

v~ 


'x\2 


...(even  function) 


(5.19) 


This  choice  is  prompted  by  the  dominant  TEM  mode  behavior  characterized  by  an 
even  longitudinal  current  and  an  odd  transverse  current  [17].  An  additional  advantage 
of  this  selection  is  that  a  closed  form  solution  for  fma(0  and  5w(0,  which  are  given 
by 

S-a  T2m+1  ( l )  \A  ~  (f)V'^G lx  if  a  =  X 

if  a  —  z 


fma(t)  = 


and 


9na  (0 


f-a  T2n+ i(f)\A  -  (f)2e-J'x5dx  if  a  —  x 


(5.20) 


f- 


a  T2"(a)  ~-3Xt.(]jX 


(5.21) 


if  a  =  z 
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can  be  found  [17].  Substitution  of(5.20)  and  (5.21)  into  (5.14)-(5.17)  leads  to 


n  r°°  ~  1 

A™  =  J  (k2  -  e)Gxx(< e,  C;  yo\yo)[J2m+M)  +  ^ (<W«)  +  J^-iK))] 

[J2n+l{aC)  +  ^  (^2n+3(a0  +  ^2n-l(a0)]^  (5.22) 

r°°  ~  1 

2^"  =  — i(  — l)ny  ^C^z(^)C;yo|2/o)[^2m+l(aO  +  +  «/2m-l(o0)] 

</2n«R  (5-23) 

45T  =  73^  J™ 

[J2n+l(a£,)  +  ^  ( "^2n+3  (®£ )  +  ^2n-l  («£)  )]  (5-24) 

/oo  ~ 

(A;2  -  C 2)GZZ{^  C;  yoll/o)  J2m(ati)J2n(am,  (5.25) 

-oo 

where  Jn  is  a  Bessel  Function  of  the  nth  order.  Since  the  integrals  in  (5.22)-(5.25) 
involve  infinite  limits,  it  is  important  to  understand  the  asymptotic  behavior  of  these 
integrands  to  address  the  issue  of  convergence.  Using  the  large  argument  approxima¬ 
tion  of  the  Bessel  function  in  conjunction  with  the  large  argument  behavior  of  the 
Green’s  function,  given  by 


Gaa  ~  -  V  a  =  x,  z  £ 


oo, 


the  asymptotic  behavior  of  the  integrands  are  found  to  be 

I(AZ)  ~  0 

I(A™)  ~  1 

I(AT")  ~  ) 

I(AT")  ~  4 


(5.26) 


(5.27) 

(5.28) 

(5.29) 

(5.30) 
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where  I  (A™*)  is  the  integrand  of  the  integral  corresponding  to  A™*.  Note,  the  desired 
computational  efficiency  cannot  be  achieved  due  to  the  slow  convergence  (less  than 
^),  thus  a  more  suitable  expansion/test  set  will  be  studied  in  the  next  section. 


5.2.3  Galerkin  Method  Using  Chehsyshev  Polynomials  of  the  First/Second  Kind. 
Improved  convergence  can  be  achieved  by  selecting  a  more  suitable  expansion/test¬ 
ing  function  set.  The  following  expansion  and  test  functions  were  selected  [17]: 


X  X 

tnx  =  enx  =  u< 2n+i)(-)Wl  -  (-)2...(odd  function)  (5.31) 

a  V  a 


T  ( — ) 

tnz  =  enz  =  -7====... (even  function)  (5.32) 

vH!F 

to  improve  computational  efficiency  while  maintaining  the  ability  to  model  current 
behavior.  New  values  of  fmz  and  gnz,  must  be  calculated  and  are  discussed  next. 


5.2.3. 1  Evaluation  of  fmz  and  gnz.  Applying  equation  (5.32)  to  the 
definition  of  fmz  and  gnz  yields 


fmZ(0  = 


f  T(2m)  (  a  ) 

f  m.)(~) 

1-aV^W- 


-Jxidx 


cos  [xff]  dx  +  j 


f  (2m)(fl)=  sjn  (x^)dx  . ..c.o.v  x  —  — 

'-a 


=0 


=  2 a  [  — (2m)('~)  cos (gxffjdx 

Jo 

=  2  aIefmM)  (5.33) 
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and 


9nz(f.) 


=  2  a 


r  wi) 
r  wi 1 

■'  r(2„)(i 


e~ix(dx 
COS  (xf)dx  —  j 


r  wi) 

/  — .  sin(x£ lax  .■.c.o.v  x  =  — 

l-a  V/R!F 


=0 


/o 

=  2a/genz(a£) 


cos  (axf)dx 


(5.34) 


where  (see  Appendix  A) 


//m,K)  =  ^K)  = 


(-l)n(f  )J2n(aO  if^O  71  =  0,1,2,3..., 
o  if  ^  =  0  71  =  1,2,3..., 

§  if(  =  0  n  =  0. 


(5.35) 


5. 2. 3. 2  Evaluation  of  fmx  and  gnx.  Applying  equation  (5.31)  to  the 
definition  of  fmx  and  gnx  yields 


fmxiO  = 


'  -a 
ra 


U{2rn+1](-)Jl-(-ye^dx 

a  V  a 


X  X  /  X 

U(2m+l)(-)\  1  -  (~)2  COs(xf)dx +j  /  C/(2m+l)(-) 

a  V  ft  / _ ri  ft 


1  —  (-)2  siri(x£)dx 
a 


=0 


=  2aj 


(5;)2  sin(ax^)dx 


=  2  ajl°fmx(af) 


(5.36) 
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and 


gmx(0  =  I  ui2n+1)(*)jl  -  {-fe-^dx 

a 


U(2n+l){-)\h-  -  (-)2  COs(x^)dx -j  /  f/(2n+l)(-) 


1  - 


,  X 


-)2sin  (x£)dx 
a 


=0 


=  —2aj  /  U(2n+i)(x)y/l  -  (x)2sin (ax£)dx 


(5.37) 


where  (see  Appendix  A) 


I fmx  (  )  I gnx  (  ®  0 


j(-l)^)^K)  if^O  n  =  0,l,2,3..., 
(o  if  e  =  0. 


5. 2. 3. 3  Asymptotical  Behavior  of  Integrands.  Again,  the  behavior  of 
the  integrands  must  be  investigated  prior  to  performing  numerical  integration.  Us¬ 
ing  (5.26)  and  the  large  argument  approximation  of  the  Bessel  function,  the  behavior 
of  the  integrands  as  £  — >  oo  is  given  by 


lirn  I(AXX) 

l 

- 

e 

(5.39) 

£— »00 

lirn  I(AXZ) 

; -oo 

1 

r s./  - 

e 

(5.40) 

lirn  I{AZX) 
£— >00 

i 

r^/  - 

e 

(5.41) 

lirn  I{AZZ) 

£— xx) 

i 

r s./  - 

e 

(5.42) 

Additionally,  when  analyzing  the  behavior  for  small  arguments,  only  the  worst  case 
scenario  (n,  m  =  0)  need  be  considered.  The  behavior  of  the  integrands  as  £  — »  0  is 
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found  to  be 


lim  I  (Axx) 
^0 

f2 

~  £ 

(5.43) 

lim  I{AXZ) 

(5.44) 

lim  I(AZX) 
^ 0 

~  £ 

(5.45) 

lim  I{AZZ) 

~  0 

(5.46) 

where  I(A™p)  is  the  integrand  of  the  integral  corresponding  to  A™p.  At  both  extremes 
(£  — >  cx)  and  £  —■ ►  0),  the  integrands  are  well  behaved  and  thus  suitable  for  numerical 
integration.  Substantial  acceleration  of  convergence  is  achieved  using  these  more  well- 
behaved  expansion  and  testing  functions.  Thus,  numerical  integration  can  be  used  to 
compute  the  MoM  matrix  elements. 


5. 2. 3. 4  Solving  the  System,  Ax  =  0.  Based  on  the  second  scenario  for 
the  testing  and  expansion  functions,  the  MoM  matrix  equation  for  the  natural  modes 


is  written  as 


Amn 

^xx 

Amn 
A xz 

Q* xn 

0 

Amn 
A zx 

Amn 
^ zz 

ttzn 

0 

where  (using  the  evenness  of  integrands) 


(5.47) 


Amn 

^xx 

poo  _ 

=  8a2/ 

Jo 

(5.48) 

Amn 
^ xz 

POO  „ 

=  8a2  j  /  -£CGgz{£,C,  Volvo) If xm{a£)Iegzn{a£)d£, 

Jo 

(5.49) 

Amn 
^ zx 

poo  „ 

=  -8a2 j  /  H(GXX(^  c,  yo\yo)I%m{<)I°9Xn(am 

Jo 

(5.50) 

Amn 
^ zz 

poo  „ 

=  8a2/  (k2-C2)Gzz^X,yo\yo)I}zm(a0le9Zn(am- 

Jo 

(5.51) 

This  matrix  equation  will  have  a  non-trivial  solution  only  if  the  determinant  is 
equal  to  zero  [36].  There  are  an  infinite  number  of  (  values  leading  to  det(A)  =  0 
and  it  is  these  values  that  represent  the  propagation  constant  of  the  natural  modes  of 
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the  stripline  structure.  Once  the  natural  mode  propagation  constants  are  found,  the 
natural  mode  currents  (i.e.,anx  and  anz )  can  be  easily  computed  using  an  eigenvector 
solution  technique.  In  this  research,  only  the  TEM  natural  mode  behavior  is  sought 
and  discussed  next. 

5. 2. 3. 5  Numerical  Results.  Except  for  the  evaluation  of  gna  and 
fma  {pt  =  x,  z )  all  integrations  were  performed  numerically.  Matlab’s  adaptive  Simp¬ 
son  quadrature  rule  was  used  exclusively  to  compute  the  matrix  elements.  A  second 
routine,  a  secant  method  root  search,  was  used  to  find  the  value  of  (  that  would  make 
the  MoM  impedance  matrix  singular.  The  null  space  of  the  system  was  found  us¬ 
ing  Matlab;  however,  numerical  errors  produced  only  a  ’’nearly  singular”  impedance 
matrix.  To  ensure  the  nullspace  was  not  empty,  an  LU  matrix  factorization  was  per¬ 
formed  and  the  value  of  the  Nth  row  and  Nth  column  of  the  resulting  reduced  matrix 
was  reset  to  zero.  To  verify  the  correct  answer  was  used  to  compute  the  current, 
the  eigenvector  corresponding  to  the  smallest  eigenvalue  was  also  used.  This  process 
provided  a  method  to  check  the  accuracy  of  the  integration  results.  As  the  precision 
of  integration  improved,  the  magnitude  of  the  smallest  eigenvalue  converged  to  zero 
as  expected.  The  numerical  results  and  effort  to  optimize  the  code  will  be  discussed 
in  the  following  sections. 

The  algorithms  designed  to  integrate  and  populate  the  impedance  matrix  can  be 
found  in  Appendix  B.  It  was  designed  to  be  accurate  and  efficient.  A  special  function 
was  used  to  select  the  upper  limit  of  integration  to  ensure  sufficient  convergence  was 
achieved.  As  seen  in  Figure  5.2,  the  the  value  for  A is  found  to  the  thousands  dec¬ 
imal  place.  At  these  large  values  of  £,  the  oscillating  integrand  is  nearly  symmetrical 
about  the  £-axis  and  the  contribution  of  the  integrand  beyond  the  limits  displayed 
in  Figure  5.2  can  be  neglected.  Similar  convergence  is  expected  for  the  remaining 
integrals,  since  they  all  converge  at  the  same  rate  (see  equations  (5.39)-(5.42)).  Effi¬ 
ciency  can  be  improved  by  exploiting  the  symmetries  of  the  impedance  matrix.  The 
and  quadrants  are  both  diagonally  symmetry.  Only  the  diagonal  and  upper 
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Figure  5.1:  Cross  Sectional  Currents  Found  on  a  Perfectly 

Conducting  Center  Conductor.  These  values  were  found  for  at 
a  frequency  of  10  MHz  for  a  stripline  with  3.46  cm  ground  plane 
spacing  and  a  center  conductor  5  cm  wide. 

off-diagonal  terms  need  to  be  calculated.  The  A ™n  and  A™n  quadrants  are  complex 
transposes;  thus,  only  values  for  one  of  these  quadrants  needs  to  be  computed.  Ex¬ 
ploiting  these  properties  reduces  the  numbers  integrations  needed  to  fill  the  impedance 
matrix  from  4(1V)4  to  2 (N)2  +  (N)  and  reduces  the  computational  time  for  hlling  the 
matrix  by  26  percent.  Furthermore,  It  will  be  shown  in  the  next  chapter  that  no 
more  than  three  expansion  terms  are  needed  to  accurately  compute  the  characteristic 
impedance,  producing  a  6x6  impedance  matrix. 

Using  the  algorithms  included  in  Appendix  B,  the  x  and  z  components  of  the 
current  and  guided  propagation  constant  were  computed.  The  currents  are  displayed 
in  Figure  5.1.  As  anticipated,  the  TEM  mode  is  being  maintained  exclusively  by 
the  longitudinal  component  of  the  current  ( Jz ).  Furthermore,  the  distribution  of  the 
current  on  the  strip  adheres  to  well  established  current  behavior  [17].  The  value  of  the 
computed  wave  number  (£  =  k  =  u^ep),  agrees  with  physical  intuition,  namely  that 
a  forward  propagating  ^-directed  wave  behaves  like  where  (  =  (3  —  ja.  For 
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2.019  2.0192  2.0194  2.0196  2.0198  2.02 

Upper  Limit  x 

Figure  5.2:  Convergence  achieved  through  special  algorithm 

(see  Appendix  B)  that  selected  the  optimal  upper  limit. 

PEC  boundaries,  there  are  no  losses  (a  =  0)  and  we  have  e^~^z\  a  pure  propagating 
wave  with  no  decay.  The  computed  currents  and  wave  number  will  be  used  in  the 
next  chapter  to  calculate  the  characteristic  impedance  of  the  lossless  stripline. 

5.3  MoM  Solution  of  the  Integral  Equation  for  Specialized  TEM  Case 

In  Chapter  4,  a  specialized  integral  equation  was  formulated  for  the  TEM  mode 
of  the  stripline  having  a  PEC  center  conductor.  It  was  shown  that  the  generalized 
EFIE  formulation  reduces  to  the  single  integral  equation  (5.52), 

a  oo 

J  J  k£Gzz(£,  k;y0  \  y0)e^xe~^x' Jsz(x')d^dx' ,  (5.52) 

—a  — oo 

under  the  condition  that  (  =  k  and  Jsx  =  0  (which  is  consistent  with  the  previous 
section).  This  integral  equation  can  be  solved  using  the  expansion  and  testing  func¬ 
tions  of  the  previous  section,  namely  enz{x)  and  tmx{x).  Performing  the  necessary 
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Figure  5.3:  Cross  Sectional  Current,  Jsz ,  obtained  using  the 
TEM  specialization  IE  for  N=3,  frequncy=10  MHz. 


expansion  and  testing  leads  to  the  matrix  equation 

N—l 

^2anzA™n  =  0...m  =  0,...,N  -1  (5.53) 

n— 0 


where 

OO 

Kz  =  f  C,  Vo  I  y0)J2n(a^J2m+2(aOdC  (5.54) 

o 

Note,  as  in  the  previous  section,  evenness  of  the  integrand  reduces  the  limits  of  inte¬ 
gration  from  (— oo,  oo)  to  (0,  oo)  and  the  leading  constants  were  removed  due  to  the 
homogeneous  nature  of  the  matrix  equation. 

It  is  shown  in  Figure  5.4  that  the  coupled  and  specialized  equation  integral 
equations  for  the  TEM  mode  current  Jsz  produce  identical  results  as  anticipated. 
However,  it  is  noted  that  the  resulting  specialized  integral  equation  MoM  matrix 
contains  N 2  elements  as  opposed  to  4iV2  elements  in  the  coupled  IE  formulation. 
This  results  a  tenfold  improvement  in  computational  efficiency. 
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Figure  5.4:  Comparison  of  both  MoM  solutions. 

5.4  Summary 

In  this  chapter,  a  MoM  solution  was  sought  for  the  IEs  formulated  in  Chapter 
4.  For  the  CIE  formulated  for  the  stripline  with  a  PEC  center  conductor,  a  Galerkin’s 
based  MoM  solution  using  Chebsyshev  polynomials  of  the  first  and  second  kind  yielded 
the  unknown  currents  and  propagation  constant.  This  approach  was  also  used  to  find 
a  very  computationally  efficient  solution  to  the  specialized  TEM  integral  equation.  In 
the  next  chapter,  the  currents  found  will  used  to  calculate  the  characteristic  impedance 
of  the  stripline. 
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VI.  Characteristic  Impedance 


The  culmination  of  the  design  process  is  to  ensure  that  the  specified  stripline  de¬ 
sign  parameters  will  produce  the  required  characteristic  impedance.  In  Chapter 
2,  the  ground  plane  spacing  was  calculated  to  provide  the  necessary  operational  fre¬ 
quency  regime.  In  addition,  an  approximate  value  for  the  center  conductor  width  was 
computed.  In  this  chapter,  an  expression  to  compute  the  characteristic  impedance 
of  the  stripline  will  be  developed  using  the  MoM  solution  obtained  in  Chapter  5. 
This  expression  will  be  used  to  test  the  accuracy  of  the  proposed  center  conductor 
width.  Any  essential  modifications  to  the  design  will  be  identified.  Furthermore,  the 
physical  understanding  of  stripline  characteristic  impedance  will  be  supplemented  by 
viewing  the  stripline  as  a  pair  of  parallel  capacitors.  Using  this  simplified  model  will 
provide  an  indication  of  factors  impacting  the  characteristic  impedance  of  a  stripline. 
Finally,  the  effects  of  center  conductor  misalignment  on  characteristic  impedance  will 
be  investigated. 

6.1  Derivation  of  Full  Wave  Expression  for  the  Characteristic  Impedance 
of  a  Stripline 

In  this  section,  an  expression  for  the  characteristic  impedance  will  be  developed 
using  the  MoM  solution  obtained  in  Chapter  5.  An  expression  for  the  characteristic 
impedance  (for  a  forward  traveling  wave)  of  a  stripline  with  PEC  ground  planes  and 
an  infinitesimally  thin  PEC  center  conductor  will  be  calculated  using  the  relation 

V 

Zo  =  j,  (6.1) 

where 

V  =  —  j  E  ■  dl . .  .C  extends  from  the  center  conductor  to  the  ground  plane  (6.2) 

c 
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and 


I  =  j)  J-  MS.  (6.3) 

CS 

By  exploiting  the  infinite  range  of  the  stripline  along  the  guiding  axis,  z  (trans¬ 
forming  to  the  ((-domain) ,  using  a  current  of  the  form  J  =  zJsz(x)5(y  —  y0)  and  se¬ 
lecting  the  most  advantageous  path  of  integration  ( dl  =  ydy ,  n  —  z  and  dS  =  dxdy), 
equations  (6.3)  and  (6.2)  can  be  greatly  simplified.  The  remainder  of  this  section  will 
outline  the  derivation  of  the  expression  used  to  compute  the  characteristic  impedance. 

The  spectral  version  of  the  characteristic  impedance  is  given  by 


V 

I 


(6.4) 


where 


V  —  —  J  E  ■  dl . .  .C  extends  from  the  center  conductor  to  the  ground  plane  (6.5) 

c 

and 

J  ■  MS.  (6.6) 

Substituting  dl  =  ydy  in  the  expression  (6.5)  yields 


V  = 


(6.7) 


where 


Ev 


9  C n 

dyujeji 


-52(0  ,y,Cn). 


(6.8) 
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Using  Leibniz’s  rule  [8],  the  voltage  can  be  found  by  evaluating  Az{Q,yXn)  at  the 
limits  of  integration.  Recalling  that  Asz(x,y,(n )  is  given  by 


A9z(x,y,Cn) 


a 

I  Gzz(x,  C;  2/0  I  lh)yj'ls(x')dx' 

—a 

a  oo 

j  j  Gzz{iX]yo\yo)eTX^e~rx'^Jls(x')didx. 

—a  —  oo 


(6.9) 


As  expected,  A*(0,  h,  (n)  (the  upper  boundary)  is  identically  zero.  The  voltage  is 
merely  ~^ASZ{0,  y„,  (n). 

The  calculation  of  Az(0,  ya,  (n)  is  very  straight  forward  when  the  results  of  the 
previous  chapter  are  used.  Using  the  expansion  of  the  current  defined  in  Chapter  5, 
the  voltage  is  given  by 


TV— 1  °9 

v  =  ^Az(0,yo,(n)  =  —  Vanz  /  &zz(£Xn,yo  I  yo)iegnz(Odt 
jixey  True  ^  J 


(6.10) 


Additionally,  the  current  reduces  to 


/  =  ancizo 


(6.11) 


where  a  is  the  half-width  of  the  center  conductor  and  azo  is  the  leading  expansion 
coefficient  for  the  longitudinal  current  Jsz  [17]. 

Substituting  the  newly  found  expression  for  the  spectral  voltage  and  current 
into  (6.4)  leads  to  the  follow  expression  for  the  characteristic  impedance  of  the  lossless 
stripline  operating  in  the  TEM  mode; 


TV— 1  00 

y-  a'nz  [  cosh  2ph  -  cosh  2py°  Te 

vr  2v/^Z7«o,./  2p  sinh  2ph 

n-l)  Q 


(6.12) 
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Figure  6.1:  The  Effects  of  Center  Conductor  Width  on  Char¬ 
acteristic  Impedance. 

This  expression  will  be  used  to  determine  the  width  of  the  stripline  field  applicator’s 
center  conductor  and  to  investigate  the  effects  of  center  conductor  misalignment. 

6.2  Validation  of  AFRL’s  Design 

Expression  (6.12)  was  evaluated  with  the  approximate  physical  parameters  sup¬ 
plied  by  AFRL.  Letting  h  =  .0173  and  a  =  .025,  yielded  a  characteristic  impedance 
of  49.76  ohms. 

While  this  is  close  to  the  desired  value,  50  ohms,  some  adjustments  are  war¬ 
ranted.  A  thinner  center  conductor  will  provide  a  characteristic  impedance  closer  to 
the  desired  value,  as  seen  in  Figure  6.1.  Re-evaluating  expression  (6.12)  with  center 
conductor  width  of  2.496  cm  leads  to  a  characteristic  impedance  of  50.09  ohms.  Note, 
more  precise  characteristic  impedance  can  be  computed  if  supported  by  fabrication 
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methods.  These  results,  as  will  be  demonstrated  in  the  next  section,  correspond  with 
physical  intuition. 


6.2.1  Physical  Interpretation.  To  ensure  reasonable  results  were  obtained 
in  the  previous  section,  it  is  useful  to  investigate  the  physical  phenomenon  of  char¬ 
acteristic  impedance.  In  this  section,  a  common  approximation  for  the  characteristic 
impedance  will  be  used  [5,31].  This  approximation  will  provide  excellent  results  for 
TEM  mode  propagation  and  give  insight  into  parameters  effecting  the  characteristic 
impedance  of  the  stripline  [1,34]. 

Consider  the  characteristic  impedance  of  the  stripline  given  by 


^o  = 


(6.13) 


where  L  is  the  inductance  per  unit  length  and  C  is  the  capacitance  per  unit  length  [5]. 
Substituting  the  phase  velocity,  given  by 


vp  = 


1 

Vlc 


into  the  expression  (6.13)  leads  to 


(6.14) 


Z0  = 


1 

vpC 


C 


(6.15) 


where  C  is  the  total  capacitance  of  the  stripline.  Neglecting  fringe  capacitance,  a 
stripline  can  be  considered  as  two  capacitors  in  parallel,  C  =  Cu  +  Ci,  where  the  total 
capacitance  is  the  sum  of  the  capacitance  between  the  center  strip  and  the  upper 
ground  plane  (Cu)  and  the  capacitance  between  the  center  conductor  and  the  lower 
ground  plane  (C\).  The  upper  and  lower  capacitances  can  be  calculated  according  to 
the  expression  given  by 

eA 

Cg  =  —  (3  —  u,l  (6.16) 

dg 
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where  dp  is  the  distance  between  the  the  center  conductor  and  the  corresponding 
ground  plane,  A  is  the  area  of  the  plates  and  e  is  the  permittivity  of  the  material 
between  the  center  conductor  and  the  ground  planes,  ft  is  clear  that  variations  in 
the  spacing  of  the  ground  planes  (h)  and  the  width  of  the  center  conductor  (a)  will 
influence  the  characteristic  impedance.  Increasing  the  width  of  the  center  conductor 
increases  the  area  of  the  capacitors  and  will  decrease  the  characteristic  impedance, 
as  seen  in  Figure  6.1.  Increasing  the  spacing  between  the  ground  planes  increases 
the  distance  between  the  plates  of  the  capacitors  and  the  characteristic  impedance 
increases.  Insight  gained  in  this  section  can  also  provide  insight  into  the  effects  of 
center  conductor  misalignment. 

6.3  Effects  of  Center  Conductor  Misalignment 

Building  on  insight  gained  in  the  previous  section,  the  effects  of  center  conduc¬ 
tor  misalignment  will  be  investigated.  First,  the  capacitor  model  will  be  adjusted  to 
accommodate  center  conductor  misalignment.  The  physical  insight  gained  through 
this  process  will  be  used  to  verify  results  obtained  by  computing  the  characteristic 
impedance  of  the  stripline  applicator  with  various  levels  of  center  conductor  misalign¬ 
ment.  Finally,  further  verification  will  be  accomplished  through  CEM  simulation. 

To  accommodate  center  conductor  misalignment,  the  total  capacitance  of  the 
stripline  can  be  written  as 


C 


Cu  +  c 
eA 


+ 


eA 


h-y  o  h  +  y0 

,,  1  1 
eA( 


h  —  yo  h  +  y0 
eA ,  1  1 


hKl-y0/h  l  +  y0/hJ' 


(6.17) 
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When  the  stripline  is  perfectly  aligned  (y0  =  0)  the  total  capacitance  is  given  by 


C 


2eA 

~h~ 


(6.18) 


Similarly,  for  small  degrees  of  misalignment  (y0  <<  1)  eqnation(6.19)  can  be  rewritten 
as 


C 


2eA 

“/T( 

2eA 


(yo  /hf 

fV0\2 


h  (1-(X>  +•• 


(6.19) 


which  reduces  to 


C 


2eA 

~h 


(6.20) 


Thus,  small  degrees  of  center  conductor  misalignment  (yo  ~  0)  will  not  strongly 
impact  the  characteristic  impedance.  Additionally,  as  the  displacement  of  the  center 
conductor  increases  the  total  capacitance  increases  and  consequently  the  characteristic 
impedance  is  driven  down.  This  will  be  demonstrated  in  the  following  section  as 
the  characteristic  impedance  of  the  stripline  is  computed  at  various  levels  of  center 
conductor  displacement. 


The  characteristic  impedance  was  computed  using  expression  (6.12)  evaluated  at 
various  values  of  center  conductor  displacement.  When  perfectly  alignment  (y0  =  0), 
the  characteristic  impedance  is  nearly  50  ohms.  As  predicted  above,  the  characteristic 
impedance  of  the  stripline  decreases  as  the  center  conductor  migrates  to  the  upper 
ground  plane  (see  Figures  6.2).  However,  the  characteristic  drops  more  rapidly  than 
expected,  as  seen  in  Figure  6.2.  Examination  of  the  electric  field  for  the  stripline  with 
the  center  conductor  located  at  yQ  =  .5 h  (see  Figure  6.3)  reveals  that  the  electric 
field  lines  fringe  from  the  bottom  of  the  center  conductor  to  the  upper  ground  plane, 
effectively  increasing  the  area  of  the  fringe  capacitors.  Thus,  the  parallel  capacitor 
model  is  limited  in  its  ability  to  predict  the  effects  of  center  conductor  misalignment. 
Due  to  the  discrepancies  between  the  expected  and  computed  values  of  characteristic 
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impedance,  additional  verification  of  the  effects  of  center  conductor  misalignment  was 
obtained  through  Microwave  Studio  simulation  (see  Figure  6.4)  [35].  The  simulated 
and  computed  results  agreed  (see  Figure  6.2).  Thus  with  confidence,  expression  (6.12) 
predicts  the  effects  of  center  conductor  misalignment. 

Expression  (6.12)  was  also  evaluated  using  various  numbers  of  expansion  terms. 
As  seen  in  Figure  6.5,  it  is  accurate  when  as  few  as  three  expansion  terms  are  used.  The 
results  for  the  characteristic  impedance  using  more  than  three  expansion  terms  are 
identical.  While  using  two  expansion  terms  provides  a  close  solution,  greater  accuracy 
can  be  gained  by  adding  more  terms.  Using  a  single  expansion  term  (N  =  1)  or  trying 
to  model  the  current  with  a  zeroth  order  Chebsyshev  polynomial  (constants  for  both 
the  first  and  second  kind  of  Chebsyshev  polynomials)  results  in  a  poor  solution. 

6-4  conclusion 

In  this  chapter,  a  full  wave  expression  was  developed  to  calculate  the  char¬ 
acteristic  impedance  of  a  striplinc  using  the  MoM  solution  obtained  in  the  previous 
chapter.  It  was  used  to  recommend  that  the  width  of  the  center  conductor  for  AFRIT s 
stripline  field  applicator  be  reduced  to  4.992  cm.  The  change  will  bring  the  charac¬ 
teristic  impedance  closer  to  the  desired  50  ohms.  Additionally,  the  effects  on  the 
characteristic  impedance  due  to  center  conductor  misalignment  were  examined.  As 
the  center  conductor  approaches  the  ground  plane  the  characteristic  impedance  drops. 
However,  small  variation  in  center  conductor  alignment  do  not  strongly  influence  the 
characteristic  impedance. 
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Location  along  x  [m] 


Figure  6.3:  Effects  of  Center  Conductor  Misalignment  of  Elec¬ 
tric  Field. 
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Type 

E-Field 

Mode  type 

TEH 

Accuracy 

2.742S7e-008 

Beta 

84.0566  1/m 

Imp.  ZPI 

50.0683  Ohms  (re) 

Maue  Inp. 

374.44  Ohms 

Line  Inp. 

50.0683  Ohms 

Plane  at  z 

“30 

Maximum-2d 

4377.36  u/m  at  2.50764 

/  -0.00856388  /  -30 

Frequency 

H  (change  with  up/doun 

keys) 

Phase 

B  degrees 

Figure  6.4: 
Results. 


Microwave  Studio  CAD  Model  and  Simulation 


0.002  0.004  0.006  0.008  0.01  0.012  0.014  0 


Center  Conductor  Displacement  [m] 

Figure  6.5:  The  Effects  of  Number  of  Expansion  Terms  (N) 
on  Characteristic  Impedance. 


VII.  Conclusion 


In  conclusion,  this  research  reviewed  the  approximate  design  parameters  proposed 
for  the  design  of  a  stripline  that  will  be  used  to  test  the  electromagnetic  proper¬ 
ties  of  materials  at  high  temperature  and  low  frequency.  It  began  by  summarizing 
fundamental  concepts  from  transmission  theory.  The  properties  of  stripline  were  ex¬ 
amined.  Additionally,  the  stripline  design  process  and  a  few  stripline  applications 
were  presented.  It  also  briefly  reviewed  recent  developments  in  stripline  fabrication. 
Next,  the  vector  potential  parallel-plate  Green’s  function  for  an  electric  current  was 
developed  using  the  spectral  and  spacial  vector  potential  boundary  conditions  for  the 
parallel-plate  wave  guide.  The  subsequent  forced  vector  potential  wave  equation  was 
solved  by  decomposing  the  equation  into  principle  and  reflective  parts.  Each  part 
was  then  solved  independently.  Complex  analysis  was  used  to  find  the  solution  to 
the  principle  part  and  the  reflected  part  was  readily  identified  as  having  a  harmonic 
solution.  Boundary  conditions  were  applied  to  the  total  solution  at  the  outer  conduc¬ 
tors  to  find  the  Green’s  function.  Next,  a  compact  hyperbolic  version  of  the  Green’s 
function  was  identified  and  used  to  formulate  an  integral  equation  for  a  stripline.  A 
pair  of  coupled  integral  equations  were  formulated  for  a  stripline  having  a  perfect  cen¬ 
ter  conductor.  Furthermore,  it  was  shown  that  the  coupled  integral  equation  reduce 
to  a  single  integral  equation  for  the  TEM  mode.  This  specialization  naturally  led 
to  a  more  computationally  efficient  formulation.  Next,  the  solution  of  these  integral 
equations  was  obtained  employing  a  Galerkin’s-based  MoM  solution  using  Chebsy- 
shev  polynomials  as  expansion  and  testing  functions.  The  unknown  center  conductor 
currents  and  guiding  axis  propagation  constant  were  then  used  to  formulate  an  ex¬ 
pression  for  the  characteristic  impedance.  The  expression  was  used  to  analyze  the 
approximate  design  parameters  of  AFRL’s  stripline  field  applicator  and  examine  the 
effects  of  center  conductor  misalignment.  It  was  used  to  recommend  that  the  center 
conductor  for  AFRL’s  striplinc  held  applicator  be  reduced  to  4.992  cm.  The  change 
will  bring  the  characteristic  impedance  closer  to  the  desired  50  ohms.  Additionally,  as 
the  center  conductor  approaches  the  ground  plane  the  characteristic  impedance  drops. 
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However,  small  variation  in  center  conductor  alignment  do  not  strongly  influence  the 
characteristic  impedance. 

7. 1  Contributions 

This  design  verification  and  characteristic  impedance  understanding  will  pave 
the  way  in  giving  Air  Force  Research  Laboratory  the  capability  of  measuring  shielding 
materials  at  low  frequencies  and  high  temperatures. 

7.2  Future  Research 

During  this  research,  much  was  learned  about  the  behavior  of  striplines.  In  spite 
of  this,  much  work  in  this  area  remains.  Adding  thickness  to  the  center  conductor, 
extending  the  research  to  striplines  that  are  variant  along  the  guiding  axis,  replacing 
PEC  with  lossy  high  temperature  alloys  and  finding  an  analytical  solution  for  the 
integrals  encountered  while  obtaining  the  method  of  moments  solution  would  be  very 
interesting. 

Literature  indicates  that  replacing  the  infinitesimally  thin  center  conductor  with 
an  actual  center  conductor  with  finite  thickness  will  improve  characteristic  impedance 
calculation  accuracy  [10]. 

This  research  was  limited  to  uniform  center  conductor  misalignment.  Utilizing 
non-uniform  transmission  line  theory  where  the  capacitance  is  dependent  on  position 
may  yield  insight  into  non-uniform  center  conductor  alignment. 

Facilitating  high  temperature  material  characterization  measurements  necessi¬ 
tates  use  of  high  temperature  alloys  in  stripline  construction.  The  conductivity  of  high 
temperatures  alloys  is  significantly  lower  than  the  conductivity  of  traditional  room 
temperature  conductors,  as  seen  in  Table  7.1.  The  effects  of  a  lossy  center  conductor 
on  the  characteristic  impedance  and  propagating  characteristics  of  the  stripline  should 
be  investigated  to  improve  the  accuracy  of  high  temperature  material  characterization 
measurements. 
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Table  7.1:  The  conductivity  of  high  temperatures 

alloys  is  significantly  lower  than  the  conductivity  of 
traditional  conductors. 


Metal 

Conductivity 

Copper 

5.9  x  107 

H-214*(  at  2200  deg  F) 

8.1  x  105 

H-230*(at  1800  deg  F) 

7.9  x  105 

HR-120*  (at  2200  deg  F) 

7.7  x  105 

H-282*(at  1800  deg  F) 

7.6  x  105 

*  High  temperature  alloys  produced  by  Haynes  International,  Inc. 

There  has  also  been  success  in  finding  analytical  solutions  to  integral  equations 
resulting  from  similar  geometries.  [12,13,38]  An  analytical  solution  using  a  pole  series 
expansion  or  special  function  could  dramatically  increase  computational  efficiency. 
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Appendix  A.  Chebsyshev  Polynomials 


This  appendix  lists  the  expressions  and  formulas  used  in  evaluating  integrals 
containing  Chebsyshev  polynomials  encountered  in  the  MoM  solution  of  Integral 
Equations.  The  main  reference  is  [7,17,33]. 

Throughout  the  appendix,  Tn(x)  will  specify  a  Chebsyshev  polynomials  of  the 
first  kind  of  order  n.  Likewise,  Un(x)  will  indicate  a  Chebsyshev  polynomials  of  the 
second  kind  of  order  n. 

Chebsyshev  polynomials  of  the  first  kind  were  generated  using  the  following 
recursive  relations,  given  as 


Tn(x)  =  2xTn-i(x)  -  Tn-2(x) 


(A.l) 


A  few  Chevbsyshev  polynomials  of  the  first  kind  can  be  given  as 


T0(x)  =1 

Ti(x)  =x 

T2(x)  =2x 2  —  1 

T:i(x)  =Ax3  —  3x 

T4(x)  =8x4  —  8x2  +  1 

T5(x)  =16a:5  —  20x3  +  5x 


(A.2) 


(A.3) 


Chebsyshev  polynomials  of  the  second  kind  were  generated  using  the  following 
recursive  relations,  given  as 


Un{x)  =  2 xUn-iix)  -  Un_2(x) 


(A.4) 
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A  few  Chevbsyshev  polynomials  of  the  second  kind  can  be  given  as 


U0{x)  =1 
Ui(x)  =2x 

U2(x)  =4x2  -  1  (A. 5) 

U:i(x)  =8x 3  —  4x 
U4(x)  =16x4  —  12x2  +  1 
U5(x)  =32x5  —  32x3  +  6a; 

(A.6) 


A.l  Preliminary  Expressions 

The  follow  expressions  are  used  to  evaluate  the  integrals  involving  Chebsyshev 
Polynomials  encountered  in  Chapter  5 


Tn(-x )  =  {-l)nTn(x) 

(A. 7) 

eJX  =  cos  x  +  j  sin  x 

(A. 8) 

Tn(cos9)  =  cos  n9 

(A. 9) 

sin2  9  =  1  —  cos2  9 

(A. 10) 

.  9  .  1  —  cos  29 

sir  9  = - 

2 

(A. 11) 

1 

cos  A  cos  B  =  -  [cos(A  +  B)  +  cos(  A  —  B)\ 

(A. 12) 

1 

1 

(A. 13) 

sin  a;  = - 

2  j 

n  Jn(x)  =  [  cos  n9ejxcosed9 

\J)n  Jo 

(A. 14) 

Jn(aejmn)  =  ejmn7r  Jn(a) . . .  m  =  integer 

(A. 15) 
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A. 2  Integrals  with  Chebsyshev  Polynomials 

Some  of  the  integrals  involving  Chebsyshev  Polynomials  encountered  in  MoM 
solutions.  Proofs  to  these  integrals  are  available  in  literature  [7,17,33]. 

rn(a)  if  a  f  0  and  n  —  0,  2, 4  . . . 

if  a  —  0  and  n  —  0  (A.  16) 

if  a  —  0  and  n  —  2,4... 


Tn{x) 
\/l  —  x 


cos  axdx  = 


-l)n/\n/2  )J 


<! 


l 

jun(x) sill  axdx 
0 


0 


if  a  xfx  0  and  n  =  1,  3,  5 

if  d  —  0  and  n  =  1,  3,  5 
(A. 17) 


A. 3  Asymptotic  Form  of  Integrals  Involving  Tn(x)  and  Un(x) 
The  above  integrals  have  the  following  asymptotic  behavior 


lim 

a— >  oo 


T 

-I  n. 


X 


:  COS 


X- 


[ax)dx 


cos[a  —  |  (n+  |) 


,1/2 


.  n  —  0,  2, 4 . 


(A.18) 


-...ii  =  1,3,5...  (A. 19) 

0 

A. f  Evaluation  of  Integrals  Encountered  in  Chapter  5 

The  following  integrals  were  encountered  will  trying  to  find  a  MoM  solution 
to  CIE  formulated  in  Chapter  4.  These  integrals  involve  the  Galerkin  based  MoM 
solution  using  Chebsyshev  polynomials  of  the  first  kind  as  expansion  functions. 


Jim  /  Un(x)Vl  —  x2  sin  axdx 


cos 


a  -  |  [n  +  | 
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A. 4-1  Evaluating  fmx(0  and  gnx(0  of  The  First  Expansion/Test  Function  Set. 


/„«({)=  /  r2m+1(-),/l  -  (-fe^dx 
a  V  a 


(A. 20) 


expanding  the  exponential  using  (A. 8)  yields 


fmx(0  =  f  T2m+i(-)\  l  -  (~)2cosxfdx  +  j  (  T2m+i(-)\  1-  (~)2  sin  xfdx 
J —a  a  V  a  J_a  a  V  a 

(A.21) 

The  first  term  on  the  right  hand  side  of  the  expression  above  vanishes  because  the 
integrand  is  odd  leading  to 


fmx{£)=j  /  T2m+1(-)Tl- (-fsinxfdx 
/  „  a  V  a. 


(A. 22) 


performing  a  change  of  variable,  where  x  —  -  leads  to 


fmx(£,)  =  aj  /  T2m+i{x)\/l  -  (x)2  sin  axfdx 


(A. 23) 


performing  a  second  change  of  variable  where  x  =  cos  9  yields 


fmx(0  —  ~ aJ  /  rE2.m+\ (cos  9 )  \J  1  —  (cos  0)2  sin  9  sin(o  cos  9ff)d9  (A. 24) 


which  can  be  rewritten  as 


fmxiO  =  aj  I  T2m+i  (cos  9)  \fl  —  (cos  9)2  sin  9  sin(a  cos  9f)9 
Jo 


(A. 25) 


substituting  (A.  16)  into  the  equation  above  and  using  (A.  10)  leads  to 


fmx{ 0  =  aj  /  cos[(2m  +  1)0]  sin2  9  sin(a  cos  9f)d9 
Jo 


(A. 26) 
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expanding  the  sin2  0  according  to  (A.  11)  leads  to 


fmx(0  —  —  /  cos  [(2m  +  1)0](1  —  cos  2 0)  sin  (a  cos  9£)d9 
2  Jo 

which  can  be  rewritten  as 

)s[(2m  +  1)0]  sin(acos0£)d0 
cos[(2m  +  1)6}  cos  26  sin(a  cos  6£)d,6 

using  the  identity  (A.  12)  will  yield 

*  pTV 

fmx(0  =  yf  /  cos[(2m  +  1)6}  sin(acos0£)d0 
2  Jo 

cos[(2m  +  3)0]  sin(acos  6£)d,6 
cos[(2m  —  1)0]  sin(a  cos  6£)dd 

expanding  the  sin  function  according  to  the  definition  (A.  13)  leads  to 

fmxiO  =  7  [\os[(2m+l)6]e^acos^d6 

cos[(2  m  +  l)6}e-j{acosd^d6 
cos  [(2  m  +  3  )6}ej{acose^d6 
cos  [(2m  +  3)6}e~j('acos0^d6 
cos [(2m  -  l)0]ei(acos^)00 
cos[(2m-  l)6}e-j(acose^d6 


(A. 27) 


(A. 28) 


(A. 29) 


(A. 30) 
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using  (A.  14)  leads  to 


fmAO  = 


T  (  C\ 

4? _ j^2m+l  Am+1  VflS  ) 

a7T  T  (  C\ 

4(_j)2m+lJ2m+l(‘  a^> 

a7T  T  (  C\ 

8(—j)  2m+3  J'2m+3  ' a ^ ' 

07T  , 

g7 _ j'^2m+3  "2m+3 1 — aS  ) 

a7F  7  / 

g  ^ _ j  ^  2m  —  1  “2m— 1  v  ’  ) 

OTT  ,  , 

+  g(_j)2m-l  ^2m- 1  ( ' ~a£) 1 & 


IX 


which,  using  (A.  15)  can  be  rewritten  as 


(A.31) 


J  Q/ft  \ 

fmxifC)  ^  j\ 2m  [Aim+l  (®£)  A  ~  ( '^2m+3(®0  A  i(o£))] 

a  similar  process  processes  yields  the  following  results  for  gnx{f)  where 


(A. 32) 


9nx  (0  —  2<—j\2n  [^2n+ 1(°0  +  2^2n+3^a^  ^2n-l(a0)]  (A. 33) 


A. ^.2  Evaluating  fmz(jL)  and  gnz(£)  of  The  First  Expansion/Test  Function  Set. 


fUZ)  = 


T- 


2m  v 


,ejx^dx 


(A. 34) 


expanding  the  exponential  using  (A. 8)  yields 


fmz(0  = 


-  2 ml 


cos xfdx  +  j 


-  2m  v 


sin  x£dx 


(A. 35) 


The  second  term  on  the  right  hand  side  of  the  expression  above  vanishes  because  the 
integrand  is  odd  leading  to 


fmz(0 


T, 


2m  V 


:  cos xfdx 


(A. 36) 
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performing  a  change  of  variable,  where  x  —  -  leads  to 


/mz(£)  aj 


T, 


2m  v 


cos ax^dx 


(A. 37) 


using  (A.  16)  leads  to 


/mz(0 


(_l)m  (vr/2)  J2m(o^)  if  a,  £  7^  0  and  m  =  0, 1,  2  . . . 

<  |  if  a,  £  =  0  and  m  —  0 

0  if  a,  £  =  0  and  m—  1,2... 


(A. 38) 


a  similar  process  processes  yields  Identical  following  results  for  gnz(£). 

The  following  relations  were  used  to  analyze  the  behavior  of  the  integrands  as 
£  — >  oo  and  as  £  — >  0 


lim  Jn(a) 

a— >  oo 


,  1-  7T 

n  2  2 


(A. 39) 


lim 

a— >  oo 


rn(a;) 

\/l  —  x 


:  cos (ax)dx 


cos[a  —  (n  +  |)| 


d/2 


(A. 40) 
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Appendix  B.  Matlab  Code 


This  appendix  contains  all  the  code  used  to  seek  a  method  of  moments  solution 
as  described  in  Chapter5.  It  is  based  on  code  developed  by  Captain  Milo  Hyde. 


Listing  B.l:  Method  of  Moment  Solution  for  Coupled  Integral  Equation. 

(appendix3/PropagationConstantSearchChebyChebyl2.m) 

function  [prop.const  ,  iter  ,  impedMat  ,  error  ,  Jx  ,  .  .  . 

Jz]  =  Pr opagat ionCons t ant Sear chChebyCheby 1 2 ( f req  ,  h  ,  .  .  . 
y_o,tol,  W,  n,  maxlter) 

7.  Based  on  code  developed  by  Cpt  Milo  Hyde 

°l  PropagationConstantSearchChebyCheby  uses  a  secant  search 
7,  to  find  the  propagation  constant  of  a  st  r  ipl  ine  (Fig  - 1 )  .  The 
7.  algorithm  expands  and  tests  the  currents  on  the  stripline’s 
7.  center  conductor  using  weighted  Chebyshev  polynomials  of  the 
7.  1st  kind.  The  initial  guess  of  the  algorithm  is  the  free 
/  space  wavenumber ,  kO .  The  algorithm  also  returns  the  number 
7.  of  iterations  it  took  to  reach  a  solution,  the  impedance 
7.  matrix,  the  error,  the  currents  on  the  stripline’s  center 
7.  conductor,  Jx  and  Jz  ,  and  the  characteristic 
7.  impedance,  Zc  ,  assuming  quasi -TEM,  Jx  =  0. 

7. 

7.  Fig-1 

1 

1  // // / / / / / / / / / / / / // / /PEC/ // // // / / // // // // / / // // / // // / 


7.  -  y  =  h 

7.  ( eo  ,  uo ) 

7.  I--  2a  --I 

7.  y  =  o 

7. 

7. 

7.  -  y  =  -h 

7.  ////////////////////pec////////////////////////////// 

1 

7t  freq  =  frequency  of  operation  in  HERTZ 


7.  h  =  half  seperation  of  plates 

7.  tol  =  tolerance  of  iterative  root  search 

7.  W  =  half  width  of  center  conductor  in  METERS  (a) 

7.  n  =  order  of  Chebyshev  polynomial  basis  f  unct  ions 
7.  maxlter  =  max  number  of  Newton -Raphson  iterations 

7. 

warning  off  MATLAB : m_warning_end_without_block 

el=8.8541878176e-12;  7»F/m  ,  free  space  permitivity 

ul  =  4*pi*10~-7;  7»H/m  permiability 
eul = [ e 1  ul]; 

kl  =  2*pi * f req* sqrt ( el *ul )  ; 
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46  k 


=  kl; 


iter  =  2  ; 

del_k=tol *k ; 

51  k2=k+del_k ; 

[A  B]  =  Det ImpedMatr ix (k ,  0,  h,  W,  n ,  eul,  freq,  tol,y_o); 

det A  =  A  ; 

[A  B]  =  Det ImpedMatr ix (k2 ,  0,  h,  W,  n ,  eul,  freq,  tol,y_o); 
det A2=A ; 

56  while  iter<maxlter 

p=k2-detA2*(k2-k)/ (detA2-detA) ; 
if  abs(p-k)<tol 

prop_const=p ; 
impedMat  =  B ; 

61  break ; 

end 

iter=iter+l ; 
k=k2  ; 

detA=detA2 ; 

66  k2=p ; 

[A  B]  =  DetlmpedMatrix (p ,  0,  h,  W,  n,  eul,  freq,  tol,y_o); 

det A2  =  A  ; 

end 

if  iter  ==  maxlter 

71  dispC’Max  number  of  iterations  reached’); 

end 

error  =  abs  (  det  A2  *  (  k2 -k)  /  (  det  A2  -  det  A  )  )  ; 

[Jx,  Jz]  =  Plot  Current s ( impedMat  ,  W,  n)  ; 

76 

end 

function  [A,  B]  =  Det ImpedMatr ix (k ,  del_k ,  h,  W,  n,  eul,  freq, 
.  y_°) 

7. 

81  7.  [A,  B]  =  DetlmpedMatrix (k ,  del_k ,  1,  d,  W,  n,  eul,  eu2  ,  freq 
tol ) 

7. 

7,  Det  ImpedMat  r  ix  returns  the  determinant  of  the  impedance 
7.  matrix  describing  the  currents  on  the  center  conductor  of 
7.  the  stripline.  It  also  returns  the  impedance  matrix. 

86  7. 

7«  k  =  propagation  constant  of  stripline 

7«  del_k  =  propagation  constant  step  for  forward  difference  ... 
calculat ion 

7c  h=  half  the  distance  between  parallel  plates 
7c  W  =  length  of  center  conductor  in  METERS 
91  7c  n  =  order  of  Chebyshev  polynomial  basis  f unctions 

7c  eul  =  vector  containing  permittivity  and  permeability 
7c  freq  =  frequency  of  operation  in  HERTZ 
7c  tol  =  tolerance  of  iterative  root  search 


tol . . . 
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7. 


96 

el  =  eul(l);  ul  =  eul (2) ; 
kl  =  2*pi * f req* sqrt ( el *ul ) ; 

101  k  =  k  +  del_k ; 

[A,  B,  C,  D]  =  EvaluatelntsFormMatrices (k ,  n,  h,  W,  kl ,  el,  ul  ,  .. 

tol , y_o  )  ; 

B  =  [A  B; 

C  D]  ; 

A  =  det (B)  ; 

106  end 

function  [A,  B,  C,  D]  =  EvaluatelntsFormMatrices (k ,  n,  h,  W,  kl,  . 
el ,  ul  ,  .  .  . 
tol , y_o ) 

7. 

Ill  7.  [A,  B,  C,  D]  =  EvaluatelntsFormMatrices  (k ,  n,  h,  W,  kl  ,  el,  ul  , 
tol , y_o  ) 

7. 

/  Evaluat  e Int  egr alsFr omMat r i ces  evaluates  xi  numerical 
7«  integrals  and  forms  the  submatrices. 

7. 

116  7.  k  =  propagation  constant  of  stripline 

7.  n  =  order  of  Chebyshev  polynomial  basis  f  unct  ions 
7«  h  =  half  spacing  of  the  groundplanes  in  METERS 
7c  W  =  length  of  center  conductor  in  METERS 
7c  kl  =  wavenumber  of  layer 

121  7c  el  =  permittivity 

7c  ul  =  permeability 

7c  tol  =  tolerance  of  iterative  root  search 

7c  y_o  =  location  of  center  conductor 

7. 

126 

7c 7c 7o 7c 7c  Integration  Limits  and  Tolerance  7c7o7c7o7c 

subdiv  =  400;  7c  Number  of  Subdivisions  on  Xi  Integration  ... 

Interval 

[limitsXi,  tol]  =  Findlnt egr at i onLimit s (k ,  h,  kl ,  tol,  subdiv, y_o) 


7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 1 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 


131 


for  i_test  =  0:n  -  1 

for  i_basis  =  0:n  -  1 

args  =  [k  W  i_test  i_basis  h  kl  el  ul  y_o  ]  ; 
for  m  =  1 : length ( limitsXi )  -  1 

136  Alsmall=(2*le-4~2-4*le-4~5)  -  (2* 1 e - 10 ~ 2 -4* 1 e - 10 ~ 5)  ; 

Altemp(m)  =  quadl ( OEval Int A 1  , 1 imit sXi (m)  , 1 imi t sXi (m  +  1 ) 
,  tol  ,  .  .  . 

0 , args)  ; 


Blsmall  =  (2*le-4~3)  -  (2*le-10“3)  ; 
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141  Bltemp(m)  =  quadl ( OEval Int B 1 , 1 imit sXi (m) , 1 imi t sXi (m+1 ) 

,  tol  ,  .  .  . 

0 , args)  ; 

Dltemp(m)  =  quadl ( OEval Int D 1  ,  1 imit sXi (m)  , 1 imi t sXi (m  +  1 ) 
,  tol  ,  .  .  . 

0 , args)  ; 

146  end 

A ( i_t est +1 , i _bas i s + 1 )  =  sum ( Altemp) +Alsmall ; 

B ( i_t est +1 , i _bas i s + 1 )  =  sum (Bltemp) +Blsmall ; 

D ( i_t est +1 , i _bas i s + 1 )  =  sum(Dltemp); 

end 

151  end 
C  =  B  ’  ; 
end 

function  [limitsXi,  tol]  =  Findlnt egr at i onLimi t s (k ,  h,  kl ,  tol,  .. 
subdiv , v_o ) 

156  % 

7.  [limitsXi,  tol]  =  FindlntegrationLimits  (k ,  h,  kl  ,  tol,  subdiv,.. 

y_o) 

7. 

7.  FindlntegrationLimits  finds  where  the  Green’s  function 
7.  becomes  NaN  .  It  returns  the  new  xi  spectral  integral 
161  7.  integration  limits  and  returns  an  adjusted  integration 
7.  tolerance  . 

7. 

7.  k  =  propagation  constant  of  stripline 
7.  n  =  order  of  Chebyshev  polynomial  basis  f unctions 
166  7.  h  =  half  spacing  of  the  groundplanes  in  METERS 
7.  W  =  length  of  center  conductor  in  METERS 
7.  kl  =  wavenumber  of  layer 
7.  el  =  permittivity 
7.  ul  =  permeability 

171  7.  tol  =  tolerance  of  iterative  root  search 
l  y_o  =  location  of  center  conductor 

7,  subdiv  =  number  of  subdivision  on  xi  integration  interval 

7. 


176  tol  =  tol*le-6; 
xi  =  1 : 1 e5  ; 

pi  =  sqrt(xi.~2  +  k~2  -  kl~2); 


181 

Gl=( cosh (2*pl*h) -cosh (2*pl*y_o)) ./ (pi . * s inh (2*pl*h)) ; 


indexNaNl  =  f ind ( isnan (G1 ) ) ; 

186  indexZerol  =  find(Gl==0); 

if  isempty ( indexNaNl )  &&  ~ isempty ( indexZerol ) 
xi_endl  =  indexZer o 1 ( 1 )  -  1; 
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elseif  ~ isempty ( indexNaNl )  &&  isempty ( indexZerol ) 

xi_endl  =  indexNaNl (1)  -  1; 

191  elseif  ~ isempty ( indexNaNl )  &&  ~ isempty ( indexZerol ) 
if  indexNaNl  (1)  <  indexZer o 1 ( 1 ) 

xi_endl  =  indexNaNl (1)  -  1; 

else 

if  indexZer o 1 ( 1 )  ==  1 

196  xi_endl  =  indexNaNl (1)  -  1; 

else 

xi_endl  =  indexZer o 1 ( 1 )  -  1; 

end 

end 

201  else 

xi_endl  =  xi(end); 

end 

limitsXi(:,l)  =  linspace(le-4,xi(xi_endl),subdiv); 

206 

end 

function  [Jx ,  Jz]  =  PlotCurrent s ( impedMat ,  W,  n) 

7. 

211  7,  [Jx,  Jz]  =  PlotCurrent  s  ( impedMat  ,  W,  n) 

7. 

7.  PlotCurrents  plots  the  electric  currents  on  the  stripline’s 
center 

7.  conductor,  Jx  and  Jz  .  It  also  returns  Jx  and  Jz  . 

7. 

216  7.  impedMat  =  impedance  matrix  of  stripline  environment 
%  W  =  length  of  center  conductor  in  METERS 
7t  n  =  order  of  Chebyshev  polynomial  basis  f unctions 

7. 


221  warning  off  MATLAB : divideByZero 

x  =  linspace ( - W/2 , W/2  ,  1000 )  ; 
[L,U,P]  =  lu ( impedMat ) ; 

U  (  end  ,  end  )  =  0  ; 

226  b  =  inv(P)*L*U; 

as_bs  =  null(b,  ’r’)  ; 
as  =  as_bs ( 1 : n)  ; 
bs  =  as_bs (n  +  1 : end)  ; 

231  U  =  ChebyU (n) ; 

T  =  ChebyT (n) ; 


11111  Creating  Polynomial  Representation  of  Currents  11111 
236  k  =  1; 

for  i  =  1:2: length (T ( :  ,  1 ) ) 

Jz (k ,  : )  =  bs (k) *T ( i  ,  :  )  ; 
k  =  k  +  1 ; 
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241 

246 

251 

256 

261 

266 

271 

276 

281 

286 


end 

Jz  =  sum  (  Jz  )  ; 


k  =  1; 

for  i  =  2 : 2 : length (U ( :  ,1)) 

Jx (k ,  : )  =  as (k) *U(i ,  :  )  ; 
k  =  k  +  1; 

end 

Jx  =  sum  (  Jx  )  ; 

Jx  =  polyval(Jx,2*x/W)  ,*sqrt(l-(2*x/W)  .  ~  2  )  ; 

Jz  =  polyval ( Jz , 2*x/W)  . / sqrt ( 1 - (2*x/W)  .  ~2)  ; 

7. 7. 1 1 1 1 1 1 1 1 7. 7. 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 7. 1 7. 1 1 1 1 1 1 1 1 1 1 7. 7. 7. 7. 7. 7. 1 1 1 1 
7, 7. 7. 7. 7.  Getting  Rid  of  Inf  for  Normalization  7//o7«7«7o 
Jzt  emp  =  Jz  ; 

Jztemp  (1)  =  ’  ’  ; 

Jzt  emp  (  end  )  =  ’  ’  ; 

7. 7. 7. 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 7. 7. 7. 7. 7. 7. 7. 7. 7. 1 1 1 1 1 1 1 1 7. 1 1 1 

JzMag  =  abs ( Jz ) /max ( abs ( Jztemp )) ; 

JxMag  =  abs ( Jx ) /max ( abs ( Jzt emp ) )  ; 
figure ; 

plot(2*x/W, JzMag,  ’ -b ’  ,  ’Linewidth1  ,2)  ; 
hold  on  ; 

plot(2*x/W, JxMag,  ’  - . r  ’  ,  ’Linewidth1  ,2)  ; 

xlabel ( ' Location  (2x/W)  :  , ’ f ontsize ’  ,12,  ’font weight’,  'bold'); 
ylabel (’ Relative  Surface  Current  Magnitude,  J_z  &  J_x ’,’ f ontsize 
,12,  .  .  . 

’fontweight  ’  ,  ’bold’)  ; 
legend(’J_z’  ,  ’ J_x  ’  )  ; 
if  mod(n-l,10)  ==  1 

title_str  =  [num2str (n-1)  ’~s~t  Order  Chebyshev  Polynomial  . 

Expansion ’ ] ; 
elseif  mod(n,10)  ==  2 

title_str  =  [num2str (n-1)  ’~n~d  Order  Chebyshev  Polynomial  . 

Expansion ’ ] ; 
elseif  mod(n,10)  ==  3 

title_str  =  [num2str (n-1)  ’~r~d  Order  Chebyshev  Polynomial  . 

Expansion ’ ] ; 

else 


title_str  =  [num2str (n-1)  ’~t~h  Order  Chebyshev  Polynomial 

Expansion ’ ] ; 


end 

title(title_str  ,  ’font  size’  ,12,  ’f  ontwe ight ’  ,  ’ bold ’ )  ; 
axis  tight  ; 
grid  on ; 


end 
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function  A  =  Evallnt A1 (xi ,  args) 


k  =  args  (  1)  ; 

W  =  args (2);  h  =  args (5); 
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i_test  =  args (3)  ;  i_basis  =  args  (4)  ; 
kl  =  args (6)  ; 
y_o  =  args (9)  ; 


291 

pi  =  sqrt(xi.~2  +  k~2  -  kl~2); 

Gl=( cosh (2*pl*h) -cosh (2*pl*y_o))  ./(pi  .  *  s  inh  (2*pl*h)  )  ; 

296  Il=-j*(-l)~i_basis*(((i_basis+l)*pi) ./(W*xi)) ,*(besselj (2*i_basis. 
+2 , xi *W) ) ; 

I3=j*(-l)~i_test*(((i_test+l)*pi)  ./(W*xi))  .*(besselj  (2*i_test+2,xi 
*W)  )  ; 


A  =  (kl"2  -  xi . ~2) ,*G1.*I1.*I3; 

301  end 

function  A  =  EvallntBl  (xi  ,  args) 
k  =  args  (  1)  ; 

W  =  args (2);  h  =  args (5); 

306  i_test  =  args (3)  ;  i_basis  =  args (4)  ; 
kl  =  args  (6)  ; 
y_o  =  args (9)  ; 


311  pi  =  sqrt(xi.“2  +  k~2  -  kl“2); 

Gl=( cosh (2*pl*h) -cosh (2*pl*y_o)) ./(pi . * s inh (2*pl*h) ) ; 

12  =  ( - 1 ) ~ i_bas i s *pi /2* be ssel j  (2* i_bas i s  , xi *W)  ; 

13  =  j*(-l)~i_test*(((i_test+l)*pi)  ,/(W*xi))  . *(besselj  (2*i_test+2, 
xi * W)  )  ; 

316 

A  =  -k*Gl .  *xi .*12  .*13; 
end 

function  A  =  EvallntCl (xi ,  args) 

321  k  =  args ( 1 )  ; 

W  =  args (2);  h  =  args (5); 

i_test  =  args (3)  ;  i_basis  =  args (4)  ; 

kl  =  args  (6)  ; 

y_o  =  args (9)  ; 

326 

pi  =  sqrt(xi.~2  +  k~2  -  kl~2); 

Gl=( cosh (2*pl*h) -cosh (2*pl*y_o)) ./ (pi . * s inh (2*pl*h)) ; 
Il=-j*(-l)~i_basis*(((i_basis+l)*pi)  ,/(W*xi))  . *(besselj  (2*i_basis. 
+2 , xi *W) ) ; 

I4=(-l)~i_test*pi/2*besselj  (2*i_test  ,xi*W)  ; 

331 

A  =  -k*Gl . *xi . *11 . *14 ; 
end 
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function  A  =  EvallntDl (xi ,  args) 

336  k  =  args  ( 1 )  ; 

W  =  args (2);  h  =  args (5); 

i_test  =  args (3)  ;  i_basis  =  args (4)  ; 

kl  =  args  (6)  ; 

y_o  =  args (9)  ; 

341 

pi  =  sqrt(xi.~2  +  k~2  -  kl“2); 

Gl=( cosh (2*pl*h) -cosh (2*pl*y_o) ) ,/(pl.*sinh(2*pl*h)) ; 

12  =  ( - 1 ) ~ i_bas i s *pi /2* bessel j  (2*  i_bas i s  ,  xi *W)  ; 

14  =  ( - 1 ) ~ i_t est *pi /2* bes sel j  ( 2* i_t est  , xi *W)  ; 

346 

A  =  (kl "2  -  k ~2)  * G1 . *  12 . *  14  ; 
end 
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Listing  B.2:  Method  of  Moment  Solution  for  TEM  Specialized  Case. 

(appendix3/TEMChebyChebyl2.m) 

function  [  impedMat  ,  Jz]  =  TEMChebyCheby 1 2  (freq  ,  h,  y_o,tol,  W,  n) 

7  based  on  code  developed  by  Cpt  Milo  Hyde 
3  7 
7 

7  TEMChebyCheby 1 2  finds  the  the  longitudinal  current  on  the  center 
strip 

7  of  a  stripline  (Fig-1).  It  implements  a  method  of  moment  ... 
solution  that 

7  expands  and  tests  the  currents  on  the  stripline’s  center  ... 
conductor  using 

8  7  weighted  Chebyshev  polynomials  of  the  1st  kind.  The  algorithm  . 
also 

7  returns  the  impedance  matrix ,  the  error  and  the  longitudinal  .  .  . 
current  s 

7  on  the  stripline’s  center  conductor,  Jz ;  assuming  quasi -TEM,  Jx 
=  0. 

7 

7  Fig-1 
13  7 

7  ////////////////////PEC///////////////////////////// 


7  -  y  =  h 

7  ( eo , uo ) 

7  I--  2a  --I 

18  7  y  =  0 

7 

7 

7  -  y  =  -h 

7  ////////////////////PEC////////////////////////////// 

23  7 

7  freq  =  frequency  of  operation  in  HERTZ 


7  h=half  seperation  of  plates 
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7„  tol  =  tolerance  of  iterative  root  search 
1  W  =  length  of  center  conductor  in  METERS 
28  */.  n  =  order  of  Chebyshev  polynomial  basis  f  unct  ions 
1 


warning  off  MATLAB : m_warning_end_without_block 

33  el  =  8 . 854 1878 1 76 e - 12 ;  %  F/m  (or  C2N-lm-2),  free  space  ... 
permit i vity 

ul  =  4*pi*10~-7;  °/0H/m  (or  T  m/A)  [1]  permiability 
eul=[el  ul]; 

[A  B]  =  Det ImpedMatr ix (h ,  W,  n,  y_o ,  eul  ,  freq,  tol); 

38 

impedMat  =  B; 

[Jz]  =  PlotCurrents ( impedMat ,  W,  n) ; 

end 

43 

function  [A,  B]  =  Det ImpedMatr ix (h ,  W,  n,  y_o ,  eul,  freq,  tol) 

% 

°/0  [A,  B]  =  DetlmpedMatrix  (k ,  del_k  ,  1,  d,  W,  n,  eul,  eu2  ,  freq,  . 

tol ) 

1 

48  */.  DetlmpedMatrix  returns  the  determinant  of  the  impedance  matrix 
describing 

7,  the  currents  on  the  center  conductor  of  the  stripline.  It  also 
returns 

“/,  the  impedance  matrix. 

% 

“/,  freq  =  frequency  of  operation  in  HERTZ 
53  “/,  h  =  half  seperation  of  plates 

“/,  tol  =  tolerance  of  iterative  root  search 

*/.  W  =  length  of  center  conductor  in  METERS 

*/.  n  =  order  of  Chebyshev  polynomial  basis  f  unct  ions 

% 

58 

el  =  eul(l);  ul  =  eul(2); 
k  =  2*pi * f r eq* sqrt ( e 1 *ul ) ; 

63  [B]  =  EvaluatelntegralsFormMatr ices (k ,  n,  h,  y_o,W,  el,  ul ,  tol); 

A  =  det (B)  ; 
end 

68  function  [B]  =  EvaluatelntegralsFormMatrices (k ,  n,  h,  y_o,W,  el, 

ul ,  tol ) 

•/. 

°/.  [A,  B,  C,  D]  =  EvaluatelntegralsFormMatrices  (k ,  n,  h,  y_o,W,  el 

ul ,  tol ) 

% 
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7„  Evaluat  e  Int  egr  al  sFr  omMat  r  i  ces  evaluates  xi  numerical  integrals 
and  forms 

73  */.  the  submatrices. 

% 

7.  freq  =  frequency  of  operation  in  HERTZ 
7,  h=half  seperation  of  plates 
“/,  tol  =  tolerance  of  iterative  root  search 
78  “/.  W  =  length  of  center  conductor  in  METERS 

7.  n  =  order  of  Chebyshev  polynomial  basis  f unctions 

% 


Hm  Integration  Limits  and  Tolerance  Z<,°h°/°h°h 
83  subdiv  =  400;  °/0  Number  of  Subdivisions  on  Xi  Integration  ... 

Interval 
kl=k ; 

[limitsXi,  tol]  =  Findlnt egr at i onLimit s (k ,  h,  kl ,  tol,  subdiv); 

°/.  y.  y.  y.  y.  y.  y.  y.  y.  7.  y.  7. 7.  y.  7. 7. 7.  7.  7.  7.  y.  y.  y.  7. "/.  7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7.  y.  y.  y.  y.  y.  y.  y.  y. 

°/olimitsXilarge  =  linspace  (limitsXi  (end)  ,  100000  ,  subdiv)  ; 

88  yolimitsXi  =  linspace  (le-4 , 150000  ,  subdiv)  ; 
for  i_test  =  0:n  -  1 

for  i_basis  =  0:n  -  1 

args  =  [k  W  i_test  i_basis  h  el  ul  y_o  ] ; 
for  m  =  1 : length ( limitsXi )  -  1 

93  Blsmall=(2*le-4~3)  -(2*le-10~3)  ; 

Bltemp(m)  =  quadl ( OEval Int B 1 , 1 imit sXi (m) , 1 imi t sXi (m+1 ) 


98 


tol  ,  0  ,  args  )  ; 

end 

B(i_test+1 , i_basis+l) 


sum (Bltemp) +Blsmall ; 


end 


end 
103  end 


function  [limitsXi,  tol]  =  Findlnt egr at i onLimit s (k ,  h,  kl ,  tol,  .. 
subdiv ) 

7. 

7.  [limitsXi,  tol]  =  FindlntegrationLimits  (k ,  1,  d,  kl  ,  k2  ,  el,  ul  , 

e2 ,  u2 ,  tol,  subdiv) 

108  7. 

7.  FindlntegrationLimits  finds  where  the  Green’s  function  becomes  . 
NaN.  It 

7.  returns  the  new  xi  spectral  integral  integration  limits  and  ... 
returns  an 

7.  adjusted  integration  tolerance. 

% 

113  7.  freq  =  frequency  of  operation  in  HERTZ 
“/,  h  =  half  seperation  of  plates 
%  tol  =  tolerance  of  iterative  root  search 
7.  W  =  length  of  center  conductor  in  METERS 
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118 

123 

128 

133 

138 

143 

148 

153 

158 

163 


7.  n  =  order  of  Chebyshev  polynomial  basis  f  unct  ions 

% 

7. 


tol  =  tol*le-6; 
xi  =  1 : 1 e5  ; 

pi  =  sqrt(xi.~2  +  k~2  -  kl"2); 


Gl=(cosh (2*pl*h) -cosh (2*pl*0) )  . /(pi . *  s inh (2*pl*h))  ; 


indexNaNl  =  f ind ( i snan ( G1 ) ) ; 
indexZerol  =  find(Gl==0); 

if  isempty ( indexNaNl )  &&  ~ isempty ( indexZero 1 ) 

xi_endl  =  indexZer o 1 ( 1 )  -  1; 

elseif  ~ isempty ( indexNaNl )  &&  isempty ( indexZerol ) 

xi_endl  =  indexNaNl (1)  -  1; 

elseif  ~ isempty ( indexNaNl )  &&  ~ isempty ( indexZerol ) 

if  indexNaNl (1)  <  indexZer o 1 ( 1 ) 

xi_endl  =  indexNaNl (1)  -  1; 

else 

if  indexZer o 1 ( 1 )  ==  1 

xi_endl  =  indexNaNl (1)  -  1; 

else 

xi_endl  =  indexZer o 1 ( 1 )  -  1; 

end 

end 

else 

xi_endl  =  xi(end); 

end 

limitsXi(:,l)  =  linspace(le-4,xi(xi_endl),subdiv); 
end 


function  [Jz]  =  PlotCurrents ( impedMat ,  W,  n) 

7. 

7.  [Jz]  =  PlotCurrents(impedMat,  W,  n) 

7. 

7.  PlotCurrents  plots  the  electric  currents  on  the  stripline’s 
center 

7.  conductor  ,  Jz  .  It  also  returns  Jx  and  Jz  . 

7. 

7.  impedMat  =  impedance  matrix  of  stripline  environment 

7.  W  =  length  of  center  conductor  in  METERS 

7.  n  =  order  of  Chebyshev  polynomial  basis  f unctions 

7. 


warning  off  MATLAB : divideByZero 
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168 


x  =  1 inspace ( -W/2 , W/2  ,  1000 )  ; 
[L,U,P]  =  lu ( impedMat ) ; 

U  (  end  ,  end  )  =  0  ; 
b  =  inv (P) *L*U ; 

173  as  =  null(b,’r’); 


U  =  ChebyU (n)  ; 
178  T  =  ChebyT (n) ; 


11111  Creating  Polynomial  Representation  of  Currents  11111 
k  =  1; 

183  for  i  =  1:2: length (T ( :  ,  1 ) ) 

Jz  (k  ,  :  )  =  as  (k)  *T  (  i  ,  :  )  ; 
k  =  k  +  1 ; 

end 

Jz  =  sum  (  Jz  )  ; 

188 


Jz  =  polyval(Jz,2*x/W)./sqrt(l-(2*x/W).~2); 


7. 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


193 


11111  Get t ing  Rid  of  Inf  for  Normalization  11111 
Jzt  emp  =  Jz  ; 

Jzt emp  (  1 )  =  ’  ’  ; 

Jzt  emp ( end )  =  ’  ’  ; 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


JzMag  =  abs ( Jz ) /max ( abs ( Jztemp ) )  ; 


198  set(0,  ’ def aulttextinterpreter  ’  ,  ’latex'); 

f  =  figure(2);  plot(2*x/W,  JzMag,  ’ LineWidth ’ ,  2); 

set (gca ,  ’Units’,  ’Inches’,  ’ OuterPosition ’ ,  [0  0.5  5  0.5+2]); 

set(f,  ’Units’,  ’Inches’,  ’ PaperPosition ’ ,  [115  4]); 

gr id ( ’ on  ’  )  ; 

203  xlabel (’ \ cent erl ine { Locat ion  (2x/W)}’,  ’FontSize’,  12); 

ylabel (’ Relative  Surface  Current  Magnitude,  J_z  ’,  ’FontSize’,  12)... 

i 

t it le ({ ’ \ cent erl ine { Current  through  }’;  ’ \ cent erl ine {TEM  ... 

Specialization}’},  ’FontSize’,  12); 
set (gca  ,  ’ f ontnaie ’ ,  ’times’); 

print (’ -depsc ’  ,  ’ TEMcurr ent  ’  )  ; 

208 

end 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1111  Integrals  ... 
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
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213 


1 1 1 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 1 1 1 1 1 1 1 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


function  A  =  EvallntBl (xi ,  args) 
k  =  args  (  1)  ; 

218  W  =  args (2)  ;  h  =  args (5)  ; 

i_test  =  args (3)  ;  i_basis  =  args (4)  ; 
kl  =  args  ( 1 )  ; 
y_o  =  args (8)  ; 

223 

pi  =  sqrt(xi.~2  +  k~2  -  kl~2); 

Gl  =  (  cosh  (2*pl*h) -  cosh.  (2*pl*y_o))  ./  (pi  .  *  s  inh  (2*pl*h))  ; 
I2=(-l)~i_basis*pi/2*besselj  (2*i_basis  ,xi*W)  ;7«gz 

I3=j*(-l)~i_test*(((i_test+l)*pi)  ./(W*xi))  .*(besselj  (2*i_test+2,xi... 
*  W  )  )  ;  7.  f  x 

228 

A  =  -k*Gl  .*xi .*12  .*13; 
end 


233  7. .  .  . 


7. 7. 7. 7. 1 1 1 1 1 7. 7. 7. 7. 7. 7. 7. 7. 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 7. 7. 7. 7. 7. 7. 7. 1 1 1 1 1 1 1 1 7. 1 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 1 1 1 1 1 1 1 


7. . . . 


7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 1 1 1 1 7. 7. 


7. . . . 


7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7. 7 1 1 1 1 1 


Listing  B.3:  Algorithm  for  computing  Chebsyshev  Polynomials  of  the  First  Kind 

(appendix3/ChebyT.m) 

function  [T]  =  CliebyT(n) 

11111  Generates  Chebyshev  Polynomials  of  the  1st  Kind  11111 
T  =  zeros (2*n , 2*n) ; 
for  i  =  0:2*n  -  1 
5  if  i  ==  0 

T ( 1  ,  2*n)  =  1; 
elseif  i  ==  1 

T (2 , 2*n-l)  =  1; 

else 

10  T ( i  +  1  ,  : )  =  2*circshift(T(i,:),[0  -1])  -  T ( i - 1  ,  :  )  ; 

end 


end 

7. 7. 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
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Listing  B.4:  Algorithm  for  computing  Chebsyshev  Polynomials  of  the  Second  Kind 

(appendix3/ChebyU.m) 


function  [T]  =  ChebyU(n) 

2  77777  Generates  Chebyshev  Polynomials  of  the  1st  Kind  7 7 7 7 "/o 
T  =  zeros (2*n , 2*n) ; 
for  i  =  0  :  2*n  -  1 
if  i  ==  0 

T ( 1 , 2*n)  =  1; 

7  elseif  i  ==  1 

T (2 , 2*n-l)  =  2; 

else 

T ( i  +  1  ,  : )  =  2*circshift(T(i,:),[0  -1])  -  T ( i - 1  ,  :  )  ; 

end 


12  end 

7. 7. 7. 7. 7. 7. 7  7. 7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7  7 
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